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^\ . We solve for the dynamics of a pair of spins, coupled to each other and also to an environmental 

0^ ' sea of oscillators. The environment mediates an indirect interaction between the spins, causing 

both mutual coherence effects and dissipation. This model describes a wide variety of physical 
^ ' systems, ranging from 2 coupled microscopic systems (eg., magnetic impurities, bromophores, etc), 

O I to 2 coupled macroscopic quantum systems. 

, We obtain analytic results for 3 regimes, viz., (i) The "locked" regime, where the 2 spins lock 

together; (ii) The "correlated relaxation" regime (mutually correlated incoherent relaxation); and 
(iii) The "mutual coherence" regime, with correlated damped oscillations. 
These results cover most of the parameter space of the system. 



> 

l/-^ ■ I. INTRODUCTION 

' In this paper we wish to estabhsh the dynamics of a model which describes a pair of spin-1/2 systems (ie., 2-level 
I systems), which are coupled both to each other and to an environmental sea of oscillators. Our Hamiltonian thus has 
the form 

H = H{ti,T2) + -ffosc({xk}) + X2; {Xk}) (1.1) 
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where Ti and T2 are 2 Pauli spins and the Xk, with k = l,2,...iV, are the oscillator coordinates. This model is 
interesting for many reasons, and we begin by listing some of them: 
■ (i) We shall see that the dynamics of Ti and T2 can actually be solved (once we have integrated out the environment V- 
Ch Such solvable models of quantum dissipative systems are rare. Previous examples. Laclude the "spin-boson" modeltl 
^ where a single 2-level system couples to an oscillaipr bath, the "central spin nmder'aQ (where a 2-level system couples 
to a bath of spins), quantum Brownian motionQ^El and the damped oscillatoitHl3 (where a single oscillator couples to 
an oscillator bath). Only the last two of these have a classical analogue. In both classical and quantum physics, 
an enormous amount about the physics of relaxation and decay has been learnt from such models for the study of 
'^^j ■ decoherence. However, significantly absent from these studies is the generalisation to two coupled quantum systems. 
d ' In classical physics such a generalisation introduces resonant coupling, "beats", and correlates the relaxation. We 
shall see here that the quantum generalisation has all this, as well as quantum interference between the 2 systems. 
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Such phenomena are of course well-known in the context of many examples, but the advantage of Eq. (1.1) is that 
it allows us to understand how things work when the 2 systems are strongly coupled (either to each other or to their 
environment). 



(ii) The model Eq. (1.1) is particularly interesting wheji one or both spins describe the low-energy dynamics of a 
macroscopic system. In other papers under preparationE3, we have-applied the results gained herein to understand 
the tunneling dynamics of coupled Josephson junctions in a SQUIDB and coupled nanomagnets. 



Perhaps the most fundamental way in which Eq. (Id) can be applied is to discuss the quantum measurement 
processed. We are aware of no previous work which attempts to discuss both a measured system and the measuring 
apparatus on an entirely. q|Uantum level, whilst both are coupled to their environment. Previous studies (in particular 
those by Leggett et al.El'olij) deal exclusively with one or the other of these systems. In this paoer we will already 
see some startling results relevant to this problem; we will devote an entire future paper to themllj. 



(iii) Models of the type Eq. (1.1) are used widely throughout physics, particularly in condensed matter physics. 
Some disparate examples should suffice to illustrate this. A set of 2 coupled Anderson (or Kondo) impurities are 
described by Eq. (l.l)EJ, with the oscillators repsesenting low-energy conduction electron excitations. A pair of 



tunneling defects coupled to phonons or electronsiiSI, or a pair of chromophores in a large molecule, coupled to 
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phonons or polaritonJli, is almost trivially described by (^]|), as are nuclear spins in a metal0. The problem of 
coupled 2-level systems is also crucial to the understanding of many glassesEll. On a more macroscopic level, we have 
alrea dy mentioned coupled SQUID's, nanomagnets, and the measurement problem. On a more microscopic level, 



(1.1) is an obvious generalisation of Wilson's nucleus-mesodij model to 2 nucleons. All these examples except for the 



nuclear spins involve strong coupling to the environment. There are of course many examples where this coupling 
is weak (perhaps the best known is that of a pair of 2-level atoms, coupled to the EM field). Such weak-coupling 
examples are well understood, and do not require the methods used below for th eir solution. 

In this paper, we shall be concerned principally with the application of Eq. (^]^) to macroscopic systems at low 



temperature. This will lead us to study a special case of (1.1), in which H = Hq + H^^^^, where 



1 1 1 ^ 

Ho = - 2 (Aiff + A^T^) + -K^^T^T^ + - ^ m^{±l + col^^^) (1.2) 

k=l 



= 2 E(4'^e'''-''^^f + e^'^-^^fDxk (1.3) 

k=l 

where the N oscillator bath modes {xk} are assumed for simplicity (but without real loss of generality) to be momentum 
eigenstates. The Pauli matrices ti and T2 represent the 2 different spins. The two spins are coupled through a direct 
interaction term Kzz giving ferro- or antiferromagnetic correlations, depending on whether K^z is negative or positive. 
In the absence of the coupling to the environmental sea, the spins have their levels split by "tunneling" matrix elements 
Ai and A2. We work in the basis in which | 1) and | J,), eigenstates of f^, are degenerate until split by the off-diagonal 
tunneling. 

Henceforth we will refer to Hamiltonians like (1.1) as HpiscEsiji, T2; {xk}), where "PISCES" is an abbreviation 
for "Pair of Interacting Spins Coupled to an Environmental Sea" . They are an obvious generalisation of the well-known 
spin-boson Hamiltonian, which has the form 

1 1 ^ 1 ^ 

H = --At" + ef- + _ ^ ^^(^ + ^2^^') + 2^' Yl ^^^k (1.4) 

k=l k=l 

whose dynamics was discussed in detail be Leggett et alS. 

The task of the present paper is to solve for the the 2-spin reduced density matrix /5(ti, T2; t'^, t^] t), produced by 
integrating out the bath oscillators; this describes the dynamics of the 2 spins. We emphasize here that such a study 
is an essential preliminary to any attempt to understand the dynamics of any of the strongly-coupled physical systems 
listed above; we devote the whole of the present paper to elucidating this dynamics. Although the application of these 
results is reserved mostly to later papers, we will comment on this herein, where possible. 

We may (very briefly) summarize our results by remarking that there are essentially three regimes of interest for 
the PISCES model, depending on how strong is the coupling of each spin to the bath, and also how high is the 
temperature T . These are 

(i) The "Strong Coupling" regime, where the 2 spins are sufficiently correlated so they lock together, and essentially 
behave as a single spin, with a quite different tunneling matrix element and a different coupling to the bath. 

(ii) The "Correlated Relaxation" regime, where the 2 spins tunnel incoherently, but nevertheless their motion is 
correlated. 

(iii) The "Mutual Coherence" regime, where the motion of one or both spins still shows (damped) oscillations, 
which are partially correlated. 

There is also a 4*'' solvable regime which we have ignored - this is the perturbative regime, where the 2 spins are so 
weakly coupled to the bath and to each other that their dynamics can be understood by perturbation expansion in 
the interaction. What we have done in the present paper is to calculate the 2-spin density matrix in analytic form in 
the above 3 regimes. We have concentrated on the case where the spins are coupled "Ohmically" to the bath, since 
this is usually the physically interesting case. In this Ohmic case, the results can be presented in the form of a "phase 
diagram" ; see section V. 

Organisation of the paper : In section II we discuss in more detail the physical origin of the PISCES model, in 
particular for the case where the 2 spins describe the low-energy dynamics of 2 macroscopic systems. In section III 
we discuss a simple "toy" model, where 2 spins couple by a direct interaction Jgrf f| only. This model demonstrates 
correlations between the 2 spins in the absence of any dissipative environment; it also makes clear the physical meaning 
of the various density matrix elements. In section IV we set up the formal solution for p for the full PISCES model, 
using influence functional methods. Section V is the most important in the paper; therein we give the detailed analytic 
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results for p in the various regimes described above, and discuss their physical significance. We summarise our results 
in section VI. The more complicated details of the calculations are relegated to a series of Appendices. 

Those readers whose primary interest is the results for the PISCES dynamics should go directly to section V, which 
can mostly be read independently of the rest of the paper. 

II. THE MODEL 

In this section we summarize the model, and sketch the way in which the truncated PISCES Hamiltonian can arise 
from a more microscopic model. No attempt is made at a detailed derivation, since this would require that we specify 
a particular physical system - only the basic ideas are outlined. 

A. Effective Hamiltonians and Lagrangians 

To see how a PISCES Hamiltonian can arise, we start with a very simple model, with a total Lagrangian given by 

L = Lo + L,r,t + Lb (2.1) 



r ^ f -2 2 2 

Lb = - 2^mk(a;k -t^kXk 



(2.2) 



Lo = -{Miqf + Ahql) - (Fi(gi) + ^2(92)) 



(2.3) 



E 



(9ic;,' + (Z2c;,0xk + ^(— — 2 
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(2.4) 



where it is assumed that Vi{qi) and ^2(92) each describe 1-dimensional 2-well potentials, wh ich are each symmetric 
under the interchange qa — > —qa, with minima at qi = ±(701/2, q2 — ±(702/2. Note that in (2.4), the effective couplings 



c^^ and c"^^ are in general complex and out of phase with each other. Thus, for example, if the two systems in Lq in 
( p.3[ ) represent two systems at different positions Ri and R2, a commonly used form for the cjj will be = Cke*'''^", 
where k is a wave vector and a = 1,2. The form we have chosen for the coupling between the enviKonmental modes 
{xfe} and the coordinates q^ is linear in both, as in the usual Feynman-Vernon/Caldeira-LeggettacJ scheme- we 
discuss the vali dity and generality of this choice below. 

Included in (2.4) are two separate countcrterms for each potential, quadratic in qa. The purpose is, as usual, to 
restore the bare potential V{qi,q2) — Vi{qi) + V2((72), after renormalisation by the coupling to the environmental 
bosons, to its original value so that we may treat it as the physical potential. In Eq. (2.2), we have not yet truncated 
the problem to a "PISCES" problem, so that no direct interaction terms appear. This is why the counterterm 
involves 2 separate contributions. If in (2.2) there had been a term 1^12(91, 92) involving qi a nd (7 2 together in a direct 
interaction, then the counterterm would no longer be separable either. Thus, in the form (|2.3| ) the PISCES model 
only involves indire c;t in teraction between system 1 and system 2, mediated by the oscillator bath. j.. 

The Lagrangian (2.1) is a simple generalization of the single 2-well model studied by L eggett et al.El, which they 
truncated to the spin-boson model. The manoeuvres required for the truncation of (2.1) are the same; we require 



simply that the separations Vl^^' and ^^o^^ of the 2 lowest levels in wells 1 and 2 respectively, from the higher levels in 
these wells, be much greater than any energy scale (such as T) that we are interested in. If we impose an ultraviolet 

cutoff ujc on the bath modes, such that VLq 



>(2) 



(1) 



^ Ai, A2,T, then we may truncate (2T) to the form 



H — Ho + Hi, 



Hp 



(2.5) 



Hp 



- 2^mk(x; 



2 2\ 
^kXk ) 



(2.6) 



Ho = --{A,ff 



A2rf) 



(2.7) 
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(2., 



and cl,^-'). 



This is just the PISCES model of Eq. (1.2) and (|1.3D ( in (|l.2|) and (1.3), the factors goi and go2 are absorbed into 



the coupUngs 

The sphttings Ai and A2 can be calculated from Lq using WKB-instanton methods. Although this procedure is 
quite complicated and technically interesting for a spin-boson problem in which the wells are not all degenerate, we 
will not discuss it here. We will simply assume that Ai and A2 are given quantities, in the spirit of the low-energy 
effective Hamiltonian, and that they are, as usual, exponentially small compared to the energy scale of V{qi,q2) in 



(2.3) (ie., exponentially smaller than JIq^-* and fl^^' respectively) 



1(2) 



HpiscES as written in (1.2) is clearly not the most general model of this kind. A much wider range of direct 



couplings is possible; instead of i^Tz^ff f| we could use 



rrdir _ }_ \ " 7^ 'i' 



(2.9) 



We could also use more complicated indirect couplings (ie., couplings to the bath) like j'Tq X^k '-k;?^''' "^"^k, with 
/i = x,y,z and a — 1,2. In the present paper we will stick to the diagonal coupling in (1.3) and keep only the 
direct coupling in fff|. Our reasoning is as follows. Just as for the single spin-boson problem, we expect diagonal 
couplings to dominate the low-energy dynamics of the combined system, since the spins spend almost all their time in 
a diagonal state (only a fraction A/fio of their time is in a non-diagonal state, when the relevant system is tunneling 
under the barrier). In certain cases the diagonal couplings can be zero (usually for symmetry reasons), and then one 

H.C.). Here we de al o nly with the generic case of 



k_L 



must include non-diagonal couplings like ^T^^-^J^d- 

dominant diagonal couplings; we shall see in any case that the physics generated by (1.3) is complicated enough as it 



Our reason for keeping the longitudinal direct coupling in (|2.9D is connected to our choice of diagonal couplings 
between the bath and the systems. The interaction in ( |1.3D will couple ff and f| to produce a longitudinal direct 
coupling of the form ij/fffl. In fact, one may quite generally observe that in a field-theoretical context, any direct 
interaction can be viewed as the result of a coupling between fi and f2 through the high-frequency or "fast" modes 
of some dynamic field; "fast" in this context simply means much faster than the low-energy scales of interest, so that 
the interaction may be treated as quasi-instantaneous. In our case, integratin g ou t these fast modes, those with a 
frequency greater than fio then produces the longitudinal static interaction in (1.2). 

One further omission from (1.2) is an applied bias acting o n one or both of the spins, of the form Cj-fJ, for example. 
This important extra ingredient will be discussed elsewhercESliB- our present interest is in the dynamic interaction 
between the spins, generated through the bath, and most of the present paper is devoted to unravelling its effects. 



B. The Oscillator Bath Representation 



The "high-energy" Lagrangian of {2A_) to (2^) already uses an oscillator bath representation of the environment, 
ie., it already involves a low-energy truncation of some more detailfid microscopic representation of the environment. 
We briefly recall the arguments that justify this representatioroo, and when they fail. 

One begins by assuming a large number N of degrees of freedom in the environment, with couplings Jij between 



the environmental excitations |i) 
"system" obeys V, ~ 0{N-^l'^). 



and \j) of order 1/A^. One also assumes that the coupling Vj of \j) to an external 



, J ■ - K^y^-, In the PISCES model, as well as the spin-boson model, the Vj are just the Ck. There 

are of course many canonical examples of such environments, such as '^He (normal, superfluid, solid, or quantum gas), 
^He superfluid and solid, metals, semiconductors, magnets, superconductors, simple insulators, etc. Notice that any 
delocalised degrees of freedom will formally satisfy these requirements, simply because their wave-functions will be 
normalised inside a box of size ^ N. 

If we treat the couplings Vj perturbatively, up-to 2"''-order, it is always possible to represent the environment as 
a bath of oscillators (c.f., Feynman and VernoncJ, pp. 153-159). We are then bound to end up with an effective 
Lagrangian of form 



^e// — ^0 + Lint + Li 



(2.10) 



where Lb takes the form ( |2.2| ), Lo{Q, Q) is the Lagrangian of the system and 
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L^nt = - ^(Fk(Q)Xk + Gk(Q)Xk) (2.11) 
k 



(or some form related to this by a canonical transformationnj) . It is often assumed that the excursions of the system are 
small, and then Lint is written as Lmt = — X^^kXkQ, as in the PISCES and spin-boson problems. However in many 
practical situations one obviously has to go to more general coupHngs (eg., the couphng Lint — ~ X](^qS*'* ""^q 
arising from a potential probe V{r) = ^VqC^'^''^] here 6q is an environmental boson operator, which could equally 



describe fermion particle-hole pairs via 6q — X^k '^k+q'^k), and one is led to (2.11). 

Thus the question of the validity of the oscillator bath representation of the environment boils down to the appli- 
cability of the 2nd-order perturbative treatment of the Vj. It is important to understand that even if there appears 
to be some breakdown of the perturbative expansion, it may easily be curable by a suitable canonical transformation. 
Consider, eg., the Luttinger Hquid with back scattering; although the {V,} are formally ~ 0{N~^^^), they are infrared 
divergent (in Caldeira-Leggett language, they are subohmic). Nevertheless, a canonical transformation to a quite new 
set of environmental modescj reveals that this subohmic form is an artefact of the original bosonisation- the coupling 
to the new modes is a conventional ohmic one. Another salutary example is provided by the coupling of solitons to 
their environment. 

However, an environment which cannot in general be mapped to an oscillator bath is an environment of uncoupled 
(or very weakly coupled) spins, or 2-level systems. In this case one finds that the couplings of the N individual spins 
{(Tfc} to the central quantum system are not ~ 0{N~^^'^), but independent of N. In some cases they are also very 
strong, so that instead of tlia-|epvironmental dynamics being only weakly affected by the quantum system of interest, 
they are in fact slaved to itaBB. _ 

From this very brief discussion (for a much more lengthy one see ref.EZI), we see that provided we can ignore the 
effects of some "spin bath" environment (eg., nuclear or paramagnetic spins, or local defects of some kind), the 
oscillator bath description should be a good representation of the low-energy environmental degrees of freedom. 

III. DENSITY MATRIX FOR TWO COUPLED SPINS (NO ENVIRONMENT) 

Before dealing with the full PISCES model, we summarize the results for the dynamics of a much simpler problem, 
in which the 2 spins couple to each other by a static longitudinal interaction Jofff|, and no environment is present. 
In a certain sense the PISCES model can be viewed as a dissipative version of this "toy model" . 

A. Single Spin 

First recall that a single spin in a bias field e, with Hamiltonian 

H = -^f, + ^r^ (3.1) 



and with eigenvalues 



E± = ±\E\ = ±^V^T^ (3.2) 



and eigenfunctions 

^P±=A±[{E±+e)\^)-Ao\i)], (3.3) 



1/2 



A, 



(Ao) 



has a "1-spin" density matrix 



in which the wave-function ^'(i) of the system, at time t, is expanded as ^'(i) — '^"(0<^a(^)i with (f>i 
<l>2 ^ \ i)- Assuming an initial state | t), p{t) then evolves as 



(3.4) 

(3.5) 

t) and 
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Et sin Et\ 

^^ sin' Et ^sinEt J ^ ' 

In the equivalent 2-well problem, the initial "wave-packet" | t) partially oscillates between the 2 wells- the diagonal 
elements give the occupation probability of the wells, and the off-diagonal elements describe oscillatory quantum 
interference between them, which is suppressed by the bias; when e 3> Aq, the particle stays in one well, in the 
absence of any coupling to the environment. 

B. Two Spins 

Now consider 2 coupled spins, with a Hamiltonian containing a simple longitudinal interaction : 

H = -i(Aiff -fA2f|) + JorffI (3.7) 

The coupling creates a "dynamical bias" whereby each spin tends to force the other spin into one or other of its wells. 
The 2-spin density matrix takes the form 

Pr,rr,r[r^(t) = Ar,rAt)K^4i) (3-8) 

(c.f., Eq (^)),in which we order the states as {|rir2)} = {| tT)i I Ti)i I iT)j I ii)}; the amplitudes ATiT2(i) are given 
by 

— - i [cos ^l-t + cos ^l^t] + ^^0(7:;— sin Q.^t + — — sin Q.-t 
z /A1-A2 . Ai + A2 . 

I /A1-A2 . „ , , Ai + A2 . 

All = — i ["^o^ ~ "^os rj+t] + iJo[TT— sin — — — sin ) (3.9) 
where at t = we have assumed an initial state It^Tj) = | IT); and the eigenfrequencies Q.± are given by 

^l=Jl + \{^i±^2? (3.10) 



so that when Jg — 0, the amplitudes in ( p.9D reduce to 



(0) Alt Asi 

j4tt = cos cos 

TT 9 9 



.(0) Alt . A2t 

A\,' = — zcos sm 

2 2 



.(0) . . Alt A2t 

A] J — —I sm cos 

iT 9 9 



<)=sin%^sin^ (3.11) 
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From this we obtain the probability PriT^it) = I^Tir2(i)P of finding the 2-spin system in state \T1T2) at time t 
(assuming a state | tt) when t = 0), and the interference between these states. When Jg = 0, the 2 spins oscillate 
independently, with product wave functio n |ri) |T2), and p{t) has oscillation components corresponding to the sum 



n+ - 


Jo 




h 


- 


Jo 






fi^ + 


Jo 



and difference frequencies (Ai ± A2) (see ( 3.11 )). When Jq is finite, we still have sum and difference frequencies ^± 
in p(t), but the behaviour is more complicated because the wave- function is no longer a product wave- function - the 
coupling "entangles" the 2 spins. In fact the cigcnfunctions are now 

*++ = ( '-"^^ )'^' ( 7777^^1 tt) + I Ti> + I it) + ~^ ' ~M ii) 

= ( ^^7^ )'^' ( "M tt) + 1 Ti) + 1 it) + J^' ' "M ii) 
- It ) [2^tTio)^ + ' + ' + ^^^T^' 

n+ + JoY^'^ f + A2 , ,,,, , A1 + A2 
* = I — — — 

with energies 



Ai 


+ A2 




-Jo) 


Ai 


- A2 


2(f7_ 


-Jo) 


Ai 


- A2 



Ai4 


- A2 




-Jo) 


Al^ 


- A2 


2{n_ 


-Jo) 


Ai- 


- A2 



E-ir-\- — — — r^-i 



= -E^+ = (3.13) 

and the way in which the entanglement occurs depends on the magnitude of Jq and its sign. For Jg < we 

have (FM) ferromagnetic coupling; when — Jg ^ Ai,A2, the ground state j^" ) tends toward the superposition 

2-1/2^1 _|_ I antiferromagnctic (AFM) couphng, the states | ti) and | it) are emphasized and for 

Jq Ai, A2, the ground state, which is now |5'++) tends to the superposition 2^^/^(| ti) + | it))- 

Consider now the probability P^^^^{t), for a system to start in | tt) at i = 0, and finish at time t in state 
|riT2). In Figs.(0) and (||) we plot these probabilities for both strong-coupling ( |Jg| 3> Ai, A2 ) and weak-coupling 
(|Jg| ^ Ai,A2). Consider first the strong-coupling regime. For both FM and AFM coupling, if the initial state 
is I tt)i we see that the system essentially oscillates between | tt) and | ii) states with an effective frequency 
A = AiA2/(2| Jg|). The oscillations between | ti) and | J,t) are very fast, frequency ^ Jg, but their amplitude is 
reduced by a factor A/ Jg. Notice the AFM system cannot relax to some combination of the lower-energy states | ti) 
and I it) (there is no environment), a problem which does not arise (for this initial state) in the FM case, where the 



initial state | tt) is already a superposition of the 2 low-lying eigenstates j.) and l^*--) in ( 3.12 ). This situation 

is precisely analogous to the strongly biased 2-well system, where an initial state | t) in the biased system is stuck 
at high energy. In the same way, there is little coherence between the (| tt): I ii)) manifold and the (| ti), I it)) 
manifold when | Jg| is large, i.e., Jg suppresses both coherence between singlet and triplet states, and the off-diagonal 
elements in p(t). 

In the opposite limit of we ak co upling the density matrix p{t) differs little from p^^^ (t) , and there is mixing between 



all four states. From Eq. ( |3.13 ) we see that for small |Jg|, the displacement of the energy levels ^ 0(Jg) unless 



Ai = A2, in which case the 2 central levels separate linearly with Jg. In this limit 

*++ = *++ + ^ [I tt) - I ti) - I it) + I ii)) (3.14) 

= + ^ - I Ti) - I it) + I ii)] (3.15) 

= - ^ [' - I Ti) - I it) + I ii)] (3.16) 

= - ^ aTTa^ - I Ti) - I iT) + I ii)) (3.17) 
This concludes our analysis of this coupled 2-spin system. 
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IV. DENSITY MATRIX FOR THE PISCES MODEL 



In this section we set up a formal expression for the 2-spin density matrix for the PISCES system. This is done 
by first integrating out the oscillator bath modes to produce an expression in terms of an influence functional; then 
we show how the resulting expression can bej-summed over instanton paths. Some of the derivation is an obvious 
generalization of the work from Leggett et al.tl, whereas other parts are somewhat less trivial. To help readers who 
are familiar with the spin-boson model, we have adopted a notation which coincides fairly closely with that of Leggett 
et al.u, with obvious generalizations. 



A. The Influence Functional 



In order to calculate the reduced density matrix for the system of interest, we will use the well known technique 
of integrating out the environmental modes in a path integral formalism. Thus, for some general Lagrangian of the 



form (2.1), with an associated action 

S[Q, {xk}] - So[Q] + 5e™[{xk}] + S,nt[Q, {xk}] , (4.1) 
(here we denote the environmental coordinates by {xk}), the reduced density matrix p propagates according to 

p{Qf,Q'f-t)^ J dQ,dQ[Jp{Qf,Q'f,t;Q,,Q[,0)p{Q,,Q',;0) (4.2) 

We write the propagator J as a weighted double path integral: 

pq{t)=Qs f<i'{t)=Q's 

MQf,Q'f,t;Q,,Q',,0)^ D[q] D[q']e^(^ok]-So[,])p^q^q>^ (4.3) 

a double integral over the paths q{t) and q'{t), starting at i = from the positions Qi and Q[ and ending at time 
t in Q f and respectively. F[q,q'] is the famous "influence functional", whose general properties are discussed by 
Feynman and VernoncJ. It incorporates the interaction between the paths q and q' with an environment. In the case 
of a system coupled independently to a set of independent external "environmental coordinates" , we have 

F[q,q']=l[Fk[q,q'] (4.4) 
k 

with i^k the functional of the k*'* external system (which in this case is just the k*'* oscillator). 

Here we are interested in the representation of the PISCES problem either via a 4 well representation, or as 
eigenstates of t^^^t^^) (which in the truncated approximation are equivalent). The path q of the combined 2-spin 
system is written as a double path q = ((/i,(?2), one for each spin (and likewise for the path q' as well as for 



Q — (Qij Q2))- The general form of the kernel Jp in equation (4.2) is then 

Jp = y D[q,] J D[q2] J D[q[] J Ai [gi] fe] fel] ^2 fel^^fel , ^l, 92, 9^] (4.5) 

where the boundary conditions are defined by 

D[qc] = / D[q^] (4.6) 

"'9„(0)=Q„,i 

with a ~ 1,2. The factors = iAi/2, A2[g2] ~ i^2/2 are the transition amplitude associated with a single flip 

(instanton) of the relevant spin (1 or 2). In an instanton calculation A ~ iloe"'^''', where Sd is the action of the 
path that minimises Lagrange's equations in imaginary time, and fio is the preexponential factor, of order the small 
oscillation frequency of the system. F[q, q'] — F[qi, q2, q[, q!^] is again the influence functional, produced by integrating 
out the oscillator environment. 

To determine the form of F[q, q'] we need the correlation functions of the boson bath. To illustrate the general case 
let us first consider a particular one, in which the oscillator modes {xk} can be classified by momentum quantum 
number k, and the two spins are considered to be at 2 positions Ri and R2. Then defining x(r,i) = J2k^~^^ '^^^i^) 
and R — Ri — R2, we require the correlation function 
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J„/3(R, T, T) = (a;(R„, r)x(R0, t))t (4.7) 
which is conveniently written in the matrix form 

T (Vl,,\-( Ji2(R,t^)A .-o^ 

J„MR,'^)- J2(^) ) ^^■^> 

with Jq,^(R, Lj) = J dre^*'^'^ Jc(/3(R, r); as before, the indices a, (3 — 1,2, and label the two different spins. The 
diagonal terms Ja — Jaa are given directly from HpjscES (Eq. ( pT^ ) and (p^)) as 

Muj) ^ ^y^^-^^Siio ~ LUk) (4.9) 
2 ^ mk^k 

k 

whereas the inter-spin spectral function is given by 

(1) (2)* 

Ji2(R, c^) = 5 V ^^^^^DiR, uji,)S{cj - c^k) (4.10) 

k 

where Z)(R, Wk) = X^k ^'''^-^k(w) is just the propagator for the bath modes (eg., in a phonon bath it is the usual 
phonon propagator). The diagonal spectral functions Ja(w) are of course equivalent to those used by Caldeira and 
Leggett. Ji2 is similar, but note that (a) in general the function Ji2(R, w) is not a separable function of R and lu, 
and (b) it is complex (it is a retarded correlation function). 

The retardation effects begin to be important once Wk > "Uk/IRI, where is the propagation velocity of the A:*'* 
bath mode. For distances |R| <C Vk/ujk, we can ignore retardation effects, and Ji2(R, w) will be separable. In this 
paper we will always work in this limit. Typically we will be interested in interactions mediated by either phonons or 
electrons. In the case of an electronic bath, one has typically the Ohmic form 

Jap{R,Uj)=U}TJafie-^^^' (4.11) 

?7i2(R) - (r7i'?2)'/'Ve(R) (4.12) 

where ryi and r]2 are local friction coefficients of the kind discussed by Caldeira and Leggett0, and Ve(R) has the form 
Ve(R) sin^{kfR)l{kfRf, in 3 dimensions. The upper cut-off uJc in an electron bath must satisfy uJc ^ ^^o^\ ^o'\ 
as discussed in section II. Typical electronic velocities are ~ the Fermi velocity Vf ^ lO^ms"^; this means that if, 
say, r^o , then retardation effects do not become important until i? ~ 50 fim. 

For a phonon bath, one usually has the superOhmic form 



^^(R,c^)=5„;3(R)ez5 j e-"/'^^ (4.13) 

5i2(R) ~ (ffiff2)'/'V4R) (4.14) 

where m > 3 in 3 dimensions, and V0(i?) ~ (aQ/R)^, with uq a lattice constant, and gi and §2 again being local 
coupling constants. For a phonon velocity '--^ 5 x lO'^ms"^, and flo ~ IK again, we now find that retardation effects 
become important for distances R 250nm or greater. For distances |R| much less than these limits, these forms 
can be used. Again, a more precise specification of Ja/3(R, w) can be given once we fix our physical problem. 
Returning now to the influence functional, we introduce the time correlation functions 



00 



3>q(t — s)= / duj J a (lli) sin uj{t — s) (4-15) 











r„(T-s)=/ du; J„{uj) coslo{t - s)coth{Lo/2T) (4.16) 
for the diagonal elements (a = 1,2), and 
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nOO 

$12(R,T - s)= I did Ji2(R,^) sina;(r - s) (4.17) 
Jo 



nOO 

ri2(R,r-s)= / doj Ji2{R,uj) cos uj{t - s)coth{uj/2T) 
Jo 



(4.18) 



for the ofT-diagonal elements. The general form of the influence functional is now 

F[qi,q2,q[,q'2]^Fi[qi,q[] Fafe,?^] Fu[qi , q2 , q[ , 4] , (4.19) 

with Fi and F2 being the influence functionals for systems 1 and 2 alone, and F12 the part representing their interaction 
through the bath. Each path qa{t), with a ~ 1,2, is comfiosed of a series of jumps between wells, which we assume 
to be located at ±goa/2. We can then deflne the hmctionfl: 

Ur) ^ q^^[qc.{r) - q'^{T)] (4.20) 

Xo.{t) = qo^[qa{T) + qUr)] (4.21) 
so that the single spin functionals are given by 

F„[t,t'] =exp(^^^ dtj^ ds{i<^^{T - s)Ur)Xa{s) -To^ir - s)Ur)Us))^ (4-22) 
whereas the interaction functional F12 is 

F12 = exp dt ds(z$i2(R, T - s)[ei(r)x2(s) + Ur)Xi{s)] 

- ri2(R,T-s)[a(r)6(s)+6(r)ei(s)])) (4.23) 



We see that in both the single spin functional Fa and in the interaction functional F12, there is a purely reactive 
"phase" correction and a purely dissipative damping term V . The damping operates when paths depart from each 
other, ie., when both ^i(t) and ^2(3) (or vice- versa) are non-zero . We may i nterp ret these equations diagramatically 
(Fig. (H)). For both the spin-boson and PISCES problems Eq. (|4.22|) and ( 4.23 ) massively simphfy, because ((''') 



and Xai"^) may only take one of the 3 values 0,±1; each time the spin a tunnels, both and Xa change by ±1. 
Notice that these equations give Fap solely in terms of the sum and difference variables between pairs of paths; the 
paths themselves have now disappeared from the formalism. 

B. Path Integral for the Density Matrix 

The formal calculation of t he de nsi ty ni atrix eleme nts (and hence all physical properties of the system) can be 
carried out starting from (O), (45), ( 4.22 ) and ( 4.23| ). In order to do this, we must integrate over 4 separate paths 



{qiit),q2{t),qi{t),q^{t), or, equivalently 6 (0, 6 Xi X2(i) )■ 

One way of handling the path integral for J would be to write it in terms of a single path over the 16 possible 
states of qi,q2,q'i and q2- Here we shaU use a procedure analogous to that employed by Leggett et al.El for the single 
spin-boson problem. Taking advantage of the fact that the Xa and are not independent (when Xa = 0, ^q, = ±1 
and vice- versa), we write F^p as an integral over 2 paths, one for each spin, in which we distinguish only between 
diagonal "sojourn" states (^^ = 0, but Xa — il) and off-diagonal "blip" states (^^ — il whilst Xa = 0)- This has 
the diagrammatic representation shown in Fig. (^). Each such path must contain an even number of transitions, 
with action — for each transition. The transitions occur at times tj (j — l,2,...2ni) for spin 1, and times 
Uk {k = 1, 2, ...2^2) for spin 2. We may then follow the standard manoeuvre of introducing a set of "charges" which 
label the states (blip or sojourn) for each path; calling these charges ?7i_/,Cij (for spin 1) and r]2k,C'2k (for spin 2), we 
allow them to have values ±1 according to 



Xi 



(r) = E ^iJ- [^(^ - - Oir - t2,+i)] (4.24) 
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= E Ci. [S{r - t2,-i) - 0{t - ^2,)] (4.25) 

3=0 

X2(t) =^ ?72U^(^-"2fe)-e(T-W2fe+i)] (4.26) 

712 

6(t) C2fc[^(^-"2fc-l)-^(T-7/2fe)]. (4.27) 

fc=0 

The sum over the 2 paths now becomes a sum over all possible arrangements of the charges r\aj and Cafci these 
charges hving on a set of discrete times between and t (with, however, certain restrictions on the paths, and on the 
charges, depending on which density matrix element is being computed). 

The explicit expansion of the density matrix p in terms of sums over the various charge and blip ordering con- 
figurations is rather messy; to make the paper easier to follow, the detailed path integral calculations are confined 
to a series of appendices. In Appendix A the detailed structure of the complete influence functional is given, along 
with explicit expressions for the blip-blip and blip-sojourn interactions, both for bath-mediated interactions between 
different states of the same spin, and also for interspin interactions. The results can be summarised as follows. The 
single spin functionals Fa are those derived by Leggett et al.EI; they can often be treated in the "dilute blip approxi- 
mation", in which only the "blip self-energy" terms (ie., between 2 times in the same blip), and interaction between a 
blip and the sojourn immediately preceding it are included. This works, provided the interactions with the bath are 
sufficiently strong that blips are rare; essentially it assumes that the environmental decoherence is strong enough to 
severely suppress the "off-diagonal" paths in the density matrix, in which the damping term Va_(j — s) in Eq. (4.23) 
comes in. Consequently, during the waiting time between two such off-diagonal excursions or "blips" , the system 
completely "forgets" the effects of the previous blip. 

The dilute blip approximation is much more problematic when an interaction term in the Hamiltonian is trying 
to force the density matrix into off-diagonal states. This point is rather crucial when we come to the interaction 
functional F12, since the main effect of the bath turns out to be the introduction of an extra interaction e(R)ff f| 
between the 2 spins, with e(R) given by 

,-(R) ^ r dco . (4.28) 



For Ohmic dissipation, the explicit form of e(R) is 

e(R) = ai2{R)uJc (4.29) 
where ai2(R) — 7;i2(R)<Zoi'Zo2/2vr?i. For superohmic dissipation, we get 

e(R) = 2r(s) (^-^^ ffi2(R)^c (4.30) 

In the subohmic case, s < 1 so that the integral diverges and the model becomes ill-defined. It is important to notice 
that the magnitude of this bias depends on the upper cut-off frequency of the bath - this is simply because the total 
strength of the interaction involves the high-energy (adiabatic) modes. Now, this bath-induced interaction simply adds 
to the original high frequency "direct" interaction to give the total adiabatic interaction between the two spins as 

:7(R) = X,,(R) + e(R) (4.31) 

The big problem for the theory is now that J'(R) plays the role of a "mutual bias" between the 2 spins, just as Jq did 
in our toy model of section III; it turns out that most of its effects come from interactions between blips of one system 
and sojourns of the other. Unless both systems are overdamped (or both locked together, so that blips on one cannot 
overlap with sojourns of the other), then it is clear that we cannot introduce any kind of dilute blip approximation 
on a single system, independently of the other, to handle the interspin term Fi2. Nevertheless we will show, in the 
next section (and in Appendices C and D) how it is possible to evaluate F12 by generalising this approximation. We 
now turn to the detailed behaviour of the density matrix. 
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V. DYNAMICS OF THE PISCES SYSTEM 



In this section, we obtain the dynamics of the 2-spin reduced density matrix. We exphcitly calculate the 4 diagonal 
matrix elements Pr^r^i't) = p{jiT2, t] "ft, 0), which are the probabilities for the system to end up in a state \T1T2) after 
a time t, having started in a state | IT)- We will concentrate in what follows on the most interesting Ohmic case, 
where the 2 spins are coupled Ohmically to the sea of oscillators. Since it is conventional to express the strength of 
this Ohmic coupling by a dimensionless constant a, we will henceforth use the index (3 — 1, 2 to label the 2-spins (we 
do not use (3 = 1/kT in what follows). We will use temperature and frequency units in which ks = h = 1, so that 
the thermal energy T and the tunneling splitting energies are A^. 

For this Ohmic PISCES problem, we will concentrate on those parameter regimes in which analytic results can be 
obtained for -Prir2(i)- Having achieved these, it is possible to build up a fairly complete picture of the behaviour of 
PriTal^)-, in all temperature and coupling range of interest. We can identify 4 regimes in the dynamics, as shown in 
the "phase diagram " of Fig. (^. 

These 4 regimes are, in the Ohmic case, defined as follows: 

(i) The "Locked Phase" (J » T, A^/a^): In this regime, the effective coupling J is so strong that the 2 spins 
lock together, in either the states | IT) or | J,t) depending on the sign of J^. This state is very similar to the strongly- 
coupled state described in section 3; we shall see that the combined "locked spin" oscillates between | H) and | ||) 
(for FM coupling), or between | tl) and | |t) (for AFM coupling), at a renormalised frequency Ac = AiA2/|J'|. 
However these oscillations are damped - in fact the locked spin now behaves like a single spin-boson system, with a 
new coupling ac = ai -|-Q;2±2ai2 to the oscillator bath, the + (— ) corresponding to ferromagnetic (antiferromagnetic) 
coupling between the spins. 

(ii) The "Mutual Coherence" phase {Ap/ap ^ J\ J > A^): Here, the thermal energy overcomes the mutual 
coupling; nevertheless if the dissipative couplings a^'s are sufficiently small (ap <C 1), it is possible for the energy 
scale A^/a/3 to dominate even if A^ < J . In this case, even though we are dealing with a strong coupling, and the 
bath dissipation is still important, some coherence in the motion of each spin is maintained - moreover, the small J 
causes "mutual coherence" between the two spins, ie., their damped oscillations are correlated to some extent. 

(iii) The "Correlated Relaxation" or High-T phase (T Afj/ap^J). In this regime, the bath causes each spin 
to relax incoherently; however, the relaxation of the two spins is still correlated (indeed each spin relaxes in the 
time-dependent bias generated by the other). 

(iv) Finally, and much less interesting, the "perturbative regime" [J ^ A^), in which the total coupling is so weak 
that the 2 spins relax almost independently; all correlations can be handled perturbatively, and only weakly affect 
the behaviour that one calculates from the standard spin-boson model. 

We notice that the locked and correlated relaxation phases always exist; but the mutual coherence phase is much 
more delicate, and only appears if the coupling of each spin to the bath is weak (ie., ap -^1). To get a better intuitive 
feeling for the energy scales and the different regimes, it is useful to imagine how the behaviour of the PISCES system 
changes as one moves along various paths in the phase diagram. Imagine first that we start off in the lower right 
half of Fig. (0), at a very low T (much less than any other energy). Then temperature fluctuations are irrelevent, 
even if each spin is strongly coupled to the bath; the effect of the induced interaction J is simply to lock the 2 spins 
into one. This "spin complex" also couples to the bath, but now with a quite different coupling; and it can also 
tunnel between 2 orientations, but at a much lower frequency. Typically its motion at finite T will be completely 
overdamped, and it will localize as T tends to zero. However if we now raise T, the thermal fluctuations increase in 
importance, and eventually compete with Jf] the spins start to unlock from each other. This crossover regime cannot 
be handled analytically, but further increase in T causes complete unlocking- we enter the correlated relaxation phase, 
where the thermal fluctuations are strong enough to destroy any coherent dynamics of the system, but the relaxtion 
of each is still correlated with that of the other. If T ^ even these correlations become irrelevant- the 2 spins 
relax independently. 

One may make a similar traversal, this time from correlated to locked phase, this time by starting at high T but 
weak coupling (top left in the Figure), and then increasing the coupling. Much more interesting is to go to low T and 
then increase the coupling. Two situations are then possible. One is where the coupling a to the bath is not small. 
Then there is no mutual coherence regime, and we find ourselves initially in the lower right quadrant of the Figure. 
Even with large a, one may still envisage a very small (eg., if the 2 spins are very far apart). Then each behaves 
independently (with weak perurbative interactive corrections), as a spin-boson system. Increasing the coupling (ie., 
increasing their distance apart) increases the correlations, until they eventually lock; the crossover cannot be handled 
analytically. 

On the other hand if a ^ 1, so each spin is weakly coupled to the bath, one can get a situation in which their 
mutual coupling, although weak, still exceeds the renormalised tunneling splitting A*. In this case the T-fluctuations 
are not strong enough to upset their mutual correlations, or to destroy coherent behaviour on the part of each spin. 
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even if T > provided aT < A*. The spin is simply so weakly coupled to the bath that it hardly sees the T- 
fluctuations, but J causes the oscillations of the 2 spins to correlate- we get "damped beats" . Further increase in 
T decoheres these completely, and further increase in J locks the 2 spins together. This mutual coherence phase is 
clearly going to be favoured if we have 2 spins close to each other (to make big enough), yet very weakly coupled 
to the bath. 

The difficulty of setting up a PISCES model which allows a mutual coherence phase is connected with the following 
important general feature of the results. In the problem of a single spin coupled to an oscillator bath, it is known 
that one may always get some coherent behaviour provided T is low enough, and also a < 1/2. However in the 
present case of 2 spins, coherent behaviour is much harder to find; even if a is very small, coherence is destroyed if 
we "switch on" the coupling between the 2 spins, once J reaches the very small energy scale A* /a. This is because 
the "dynamic bias" on one spin, coming from the other, removes the "resonance condition" between initial and final 
states which is necessary for coherent oscillations. It shows, in essence, that mutual interactions increase decoherence 
effects. This enhanced decoherence is particularly interesting in the context of the measurement operation, or indeed 
for any attempt to look for macroscopic quantum coherence. 

A more formal understanding of these energy scales comes out of a renormalisation group treatment, to which we 
proceed first. After this the analytic results for each regime are presented. 



A. Energy Scales From the Renormalisation Group 

Although a renormalisation group analysis cannot give us the dynamics of the PISCES model, it gives considerable 
insight into the energy scales mentionad-iabove. We shall use the "poor-man" renormalisation procedure first applied to 
the Kondo problem by Anderson et al.c3, and then later to the single spin-boson problemc^l. CardjJp devised a method 
to treat the multisite problem, which was then applied to the case of two Anderson impuritiealZI and two tunneling 
impurities in metallic. As shown in Appendix B, in the case where both spins are identical (ie., where Ai = A2 = A, 
and ai — a2 = ct) the partition function for this problem is identical to that considered by Chakravarty and HirschliZI, 
and we carry over directly their scaling equations for a set of fugacities 



dy 
dhiTc 

dyp 
d\nTc 
dyAF 
dlnTr 



y{l - a) + v{yF +yAF) 

yF{l-a/2-ai2/2) + y^ (5.1) 
yAF{l-a/2 + ai2/2) + y^ 



where the dimensionless fugacity y = A/lOc and yp and yAF are the fugacities associated with simultaneous transitions 
of the 2 spins. Since a 3> ai2, the behaviour of the fugacities is determined solely by the single spin dissipation 
coefficients. These coefficients are also renormalised; their scaling is 



da 
dlnTc 

do-12 /r, 2 , 2 , 2 \ ,2 2 



a{2y'^ +yp + y\p) - 012 (y| - y\p) 

ai2{2y^ + yl + ylp) - aivl - Vaf) (5-2) 



dluTc 

Finally, there is also a scaling of the coupling energy 

The renormalisation group analysis assumes that the dimensionless quantities y, TcT and TcJ7 are all much less 
than one. The scaling has to be stopped if one of the parameters renormalises to the strong coupling regime. The 
three important scaling parameters are thus 1/A*, a renormalised tunneling matrix element, 1/T and 1/J7. Let us 
ffist set = 0. The scaling equations arc then very similar to the case of the single spin-boson system. At T = 0, 
the scaling of y is determined by a. For a > 1, y decreases and since yp and yAF growj-as y^, this decrease cannot be 
compensated and the scaling leads to y ~ 0, the well-known localisation phenomenonld. The dynamics will then be 
determined by yp or yAF- Depending on the strongest of the two, only the ferromagnetic or antiferromagnetic states 
will be occupied. For a < 1, then the scaling can stop at some r = 1/A*, a renormalised tunneling matrix element 
determined self-consistently as 
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) 



a/(l-a) 



(5.4) 



In this case, the behaviour is determined by the single tunnchng fugacity, the two spins are stiU correlated, but behave 
nearly independently. At finite temperature, the value of the renormalised tunneling matrix element is determined 
by the ratio between the energy scales 1/A* and 1/T. If T is the largest energy scale, the scaling must be terminated 
at Tc ~ l/T and the resulting tunneling matrix element is A*(r) = Ao(r/^lo)"/^^""^- 

In the presence of an interaction, the symmetry between the levels is broken and Eq. ( ^.3[ ) indicates that this 
difference will grow with the scaling, even if the fugacity grows as well (since y ^ 1 by assumption). The time it takes 
for TcJ to grow to the strong couphng phase is roughly Tj^XjJ. If both and l/T are smaller than l/^/, 

the scaling stops before the strong coupling phase is attained and the system behaves like two single spins, although 
coupled by the interaction J . However, if J is the largest energy scale, the system rapidly gets into the strong 
coupling phase and the RG equations break down. It is clear however that the new strongly coupled phase is made 
of the two spins locked together by J . The thermal fluctuations (controlled by T) and the tendency to individual 
dynamics (determined by A^) cannot compensate this strong coupling. 



Here J is so strong that a prohibitively high energy cost is incurred if a blip on one path overlaps with a sojourn 
on the other. In Appendix B this is discussed in more detail, and we demonstrate that for a strongly FM-coupled 
system which starts in initial state | tT)i the resulting dynamics is is equivalent to that of a single spin- boson system, 
with a new tunneling matrix element Ac = AiA2/|i/|, and with this "spin complex" now coupled to the environment 
with a new Ohmic coupling ac = ai + a2 + 2ai2. The new oscillation frequency is not surprising- we have already 
seen it in our toy model in the strong coupling regime. As for the new coupling to the bath, it is simply equivalent to 
the coefficient Kp p coming from the renormalisation group analysis (see Appendix B). Physically, in the case of FM 
coupling, the system oscillates between | "ft) and | J,|) at a frequency Ac, if we ignore the environmental dissipation. 
Mixing with the high-energy states | ti) and | If) is negligible (with J ^ T ^ even thermally-excited transitions, via 
absorption of bath quanta, will be exponentially rare) , and so the system reduces to the single spin-boson problem. 

A system coupled by a strong AFM interaction will have the same effective tunneling matrix element Ac, but will 
now be coupled to the environment by a dissipation coefficient ac — ol\ + oii — 1a\i. Once again, the new coupling 
can be related to a the RG analysis, in this case, it is simply the coefficient Kaf,af of Appendix B. 

Notice that the case of an AFM-coupled system which begins in the state | IT) is quite different, since this is then 
a high energy state. Without the bath, the system is frozen in this state (c.f., section 3). The coupling to the bath 
allows relaxation to the low-energy manifold (now composed of the states | ti) and | it))- Thus this case needs to 
be treated slightly differently. Intuitively one expects that the system will first try to relax to an AFM state before 
simultaneous transitions take place (since a transition from | tt) to | ii) still takes the system to a high-energy state); 
this turns out to be correct. 

Let us begin with the FM-coupled case. The equivalence to the unbiased spin-boson system means that tlie dilute- 
blip approximation can be used, and we will simply take over the results already derived by Leggett et afl. Thus, 
defining the Laplace transform of P^^it) as 



B. Locked Phase 




(5.5) 



one finds in this locked phase 




(5.6) 



The function /(A) is defined in terms of an effective tunneling matrix element A^-^-^, viz.. 



Af / ^ [r(i 



2ac)cos(7rac)]2<i-°<=' A* 



(5.7) 



Ac (Ac/fio)" 



(5.8) 



Ac 



A1A2 
J 



(5.9) 



14 



The function /(A), in the Umits appropriate to the truncation procedure (T, Ac fioi as well as ujct ^ 1) is given 

by 

/(A, T) = AJ^ j r(l - a, + X/2nT) ^^'^^^ 

where r(a;) is the Gamma fimction and where, as noted above, ac = ai + a2 + 2ai2. The corresponding behaviour 
of Pi^it) is well known: 

(i) At T = and for < 1/2, the system exhibits damped oscillations of frequency ~ A^jf. For 1/2 < etc < 1, the 
dynamics is an exponential decay superposed on an incoherent background and for ac > 1, the system is localised, 
that is P^^{t) = 1 for all times, (ii) At finite T, this oscillatory behaviour persists as long as T < Acji^fXc- For 
T > Acj-Kac, the dynamics is overdamped, and the system relaxes exponentially to P^f = 1/2. 

The various regimes are shown in Fig. (|^) . For the present problem the system is most likely to be in the overdamped 
regime; because A'jJf is so small, extremely small values of ac and/or T will be required to see underdamped 
oscillations. In the overdamped regime, the damping rate is 



A2 



2ttT 



no 

and, provided T > Ac/irac, one has the simple result 



-1 r2 



r^(ac) 
r(2ac) 



(5.11) 



Pn{t)^l + le~^'' (5.12) 

If T < Ac/ac, numerical methods must be employed if one is interested in the exact behaviour of the system. Notice 
that when ac > 1/2, an increase in temperature increases Fc, corresponding to faster relaxation. However, if etc < 1/2, 
the relaxation rate Fc (which is now much smaller) decreases with increasing T. Thus the thermal environment plays 
a different role on either side of the '"Toulouse line" ac = 1/2; for ac > 1/2, it actively retards relaxation. In the 
limit ttc <C 1, we can establish an upper bound on the relaxation rate Fc < Ac/2tt. 

Obviously, the same discussion will go through if the coupling is antiferromagnetic but with an initial | | j) or | jf) 
configuration. 

Consider now the case of an AFM coupling with an initial | IT) state. As mentioned before, the system must first 
relax to a low energy state before any simultaneous transitions can take place. In the present case of strong coupling 
this is obvious- the energy gap between the upper "FM-ordered" doublet and the lower AFM doublet is so large that 
this is by far the fastest transition, and so we do not attempt to derive this result. The initial relaxation is then 
simply determined by the relaxation of a single spin-boson system in the presence of a static bias ^ T 

Tafm = ap{Al/J){Jinof"^ (5.13) 

where f} represent the spin with the fastest relaxation rate. Once the system has relaxed to the low-energy states 
{Ti} or {it}, the analysis for the coupled dynamics can be applied. 

It has to be said that the locked phase is the one that is most likely to be encountered in physical applications. 
For microscopic models such as the problem of 2 coupled magnetic impurities, it represents the ordered state; for 
macroscopic models such as 2 coupled nanomagnets or 2 coupled SQUID 's, it will prevail unless the mutual coupling 
is very weak (ie., they are very far apart, since usually a is not small for such systems). Many quantum measurement 
schemes actually require that a strong correlation be established between apparatus and measured system at some 
point, thereby freezing the system dynamics (even if the system is otherwise able to display coherent tunneling). 



C. Correlated Relaxation Phase 



As the coupling J is reduced, or T is increased, temperature fluctuations break up the "locking of blips" which 
allowed the solution of the locked phase. When T ^ J^, we may use instead an approximation in which blips of one 
path almost always overlap with with sojourns of the other. It is important that this approach only works outside of 
the perturbative regime - but it is the strong-coupling regime we are interested in here. 

Again, the detailed evaluation of the path integrals, etc., is given in Appendix D. There we show that one can 
ignore both interactions between blips on the same path, and between blips on the 2 separate paths. This requires 
^ A/3, and thus places us directly outside the range of a perturbative treatment. Our results therefore cannot 
be necessarily extended to the limit J 0, especially if one is interested in the coherent properties of the coupled 
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system in this limit. However, provided that the two conditions T J and J A^, then our results are valid for 
any values of the dissipation coefficients a^. 

In fact one finds that the same path summation technique also works for the mutual coherence phase, and so we 
first give a general expression for the Laplace transform of -Pti(i) (again we assume an initial state | tt))- As shown 
in Appendix D, P||(A) is given in terms of the 2 functions 



50(A) = A2 / dte-^*-'?^ (*) cos\Qf\t)\ cos[Jt] 



/i0(A) 



2(//3(A + *J7) + /MA-*^)) 



/o 

tauTra 

2i 



dt( 



-\t- 



>it}. 



sin[Q^''^(i)] sm[Jt] 



(5.14) 



(5.15) 



where fp{X) is equivalent to Eq.(5.10) provided that the appropriate changes of parameters are made. We then have 



51 1 + - 
.91 



^2 

.92 1 + — 
.92 



1 



.91 ( 1 + — ) + f 1 + — 
91 J V .92 



1 

A A + gi + 52 
1 



A2 + A(5i +52) +.9152(1 -ftif) 



(5.16) 



T5i52 



1 - 



5132 



1 



A A2 + A(.gi +.92) +.91.92(1 -ft) 



Analogous expressions for the other probabilities appear in Appendix D. The dynamics is controlled by the poles of 
^^^^.^(A). There are two types of denominator; the "global denominator". 



A + 51(A) + 52(A) 

corresponds to the overall system dynamics, whereas the "individual spin" denominator, 

(A + 5i(A))(A + 52(A)) -/ii(A)/i2(A) 



(5.17) 



(5.18) 



corresponds to the individual spin dynamics. Both spins are correlated through the terms in 5/3//1/3, which essentially 
behave as tanhj7'/2T at long times. This correlation is present both in the denominators (which give the decay 
rates and oscillation frequency) and in their coefficients (controlling the weight of each term) which results in a very 
complex behaviour. 

On the other hand the limit of XPrj^r2{X) as A ^ yields the universal long time probabilities 



PTT(i ^ 00) = Puit ^ 00) = -[1 - tanh^/2T] 



(5.19) 



Pnit ^ (x.) = Pxt(^ ^ ^) = ^[l + tanh J/2r] 



(5.20) 



These results are of course nothing but the thermodynamic equilibrium expectation values, once all dynamic relaxation 
to equilibrium has ceased. 

Let us now deal with the regime T ^ J', and also assume that T > A^/a^j; this is the overdamped correlated 
relaxation regime (cf. Fig. (|)). Formally this regime is very easy to dealjwith, since we may set A = in //5(A), and 
the explicit dependence on hi and /12 in PriTiiX) can be eliminated usingH hp{Q)/gp{Q) — tanh(J'/2r). The Laplace 
inversion results in a sum of decaying exponentials for each ^^^^-^(t): 



Pn(<) = A_ + A+e-(^i+^^)* + B+e 



R+t 



(5.21) 



PuW = ^+[l-e-(^^+^^'*] + C+[e 



-R+t -R-U 



(5.22) 
(5.23) 
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P^^{t)=A+[\-e 
with the different constants defined as 



-(ri+r2)ti 



C+[e 



-R+t 



-R-t-i 



(5.24) 



A, 



-[l±tanh(J^/2T)] 



(5.25) 



^- = 4 



1 ± 



ri + r2 



tanh(J^/2T) 



(5.26) 



i?± = i(ri + r2)±i. 



(5.27) 



i?i2(T) = r2 



2rir2 [1 - 2tanh2(J/2r)] 



(5.28) 



(5.29) 



The relaxation rates Fi and T2 are [pimply the functions .gi(0) and 52(0), the relaxation rates of a spin-boson system 
in tlic presence of a bias. These are 



2^0 



2'kT' 



cosh{J/2T) 
r(2a/5) 



T{ap + ij /2TrT)\^' 



(5.30) 



Eg. ( ^.13 ) corresponds to the limit J' ^ T of this equation. In the other limiting case, J and for ^ 0(1), 



A2 r 
2Vln 



2nT 



-1 2a,3-l 



r{2ap) 



+ 0{JlTf {JIT < 1) 



However, if a/3 ^ 1 as well, then it is 



21^0 



2'kT' 



no 



aj + {J/2nTY 



{J/T,ap^ 1) 



(5.31) 



(5.32) 



Thus one can get different behaviour depending on the ratio between aj and J /T 

The relaxation times appearing in ^Vir2(0 ^ J can be quite different, and their contributions will in general 
appear in different proportions. Fig.(||) illustrates this by plotting the 4 probabilities. 

Let us now look more closely at the structure of ^'tit2(0- If = 0, so that the two system are completely 
uncorrelated, it simply decomposes into a product of the individual spin probabilities PTiT2it) = -fVi (i)^ri (0' with 
A± — B± — C-i- = 1/4, _R+ = Fi and i?_ = F2. Turning on the interaction J causes the individual probabilities 
to merge together to produce a term corresponding to the relaxation of the system as a whole (with relaxation rate 
Fi + F2, and weight A±), mixed in with correlated individual relaxation (of rates R±, and weight B± or C+). The 
mixing of the relaxation rate is independent of the sign of J] as J increase, i?+ tends to Fi + F2 while i?_ — )■ 0, 
indicating a reduction of the individual dynamics of the spins. However, the coefficients A±, B± and C+ do depend 
on the sign of the interaction, as we would expect - we have started the system in the state | It). Notice that 
the correlations of the two spins is only weakly felt in the relaxation rates R±; it represents a correction of order 
J/T << 1. The relaxation at long time will reflect this correlation more strongly, mostly through the terms A±. 

It is important to remember that all these expressions assume J <^ T (and can thus be expanded in powers of 
J /T). Using them ior J ^ T would be incorrect- one is then in the crossover region between the correlated and the 
locked phase. In this crossover, we cannot solve the problem analytically - the assumption of non-overlapping blips is 
breaking down, but they are still not fully correlated (as in the locked phase) . Physically, when J ^ T the spins are 
trying to lock, but thermal fluctuations are preventing them. 
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D. Mutual Coherence Phase 



U ap > 0(1), then we have exhausted all possible phases of the PISCES system; either the 2 spins are locked, or 
they exhibit correlated overdamped relaxation, depending on the ratio of J7 to T. 

However if ^ 1, a new possibility emerges; for temperatures such that <C T <C Ap/ajS, the thermal 
fluctuations are strong enough to prevent the locking of the blips but still not capable of destroying coherence. As 
a result, the dynamics is composed of underdamped oscillations. In the PISCES problem the oscillations of the 2 
systems will be correlated, and we can even expect decaying "beats" , between the 2 spins. 

Because we are in the non-perturbative limit J7 ^ A (even though <C 1), the oscillation frequencies are roughly 
equal to J7, with corrections of order 0{{Ap/J)'^) and O{Al{a0T)'^ /J"^). The re are 2 characteristic relaxation rates 
in the system. There is a "global" decay at a rate Tgiob ~ ai3T{Ai3/ J)"^ , and a much faster damping of the mutual 
oscillations, at a rate Tosc ~ a^T ^ Tgiob - We give more detailed expressions below. 

As noted above, expressions like ( |3.16| ) for PT1T2W a-re still valid for this regime, but an expansion of fp{X,T) in 
(5.10) (and hence of g/3(A) and /i/3(A)) in powers of A is invalid - instead we must expand in (a^ + X/2ttT) <C 1. This 
expansion gives 



A *2 

^13 



27rT 



2Q£ 



1 



r(l - 2a/3) 27rTap + A 



(5.33) 



In general the denominators in PriT2 a-re of polynomial degree 3 or 6, which renders an analytic inverse Laplace 
transformation hopeless. We can however study 2 cases. One is the case where both spins are equivalent, ie., 
Ai = A2 = A and ai = a2 = a; this is the Equivalent Spin case. We can also examine what happens if one spin is 
overdamped, and the other is underdamped; this is the Overdamped plus Underdamped case. We re-emphasize that 
in real physical situations it will be quite difficult to see mutual coherence. One possibility for macroscopic systems 
involves 2 conducting nanomagnetic particles imbedded in a semiconducting substrate; the low electronic density so 
reduces a and the mutual RKKY interactions that the equivalent spin phase is a real possibility. The "overdamped 
plus underdamped" system is interesting in view of the possibility of coupling an overdamped measuring apparatus 
to an underdamped quantum system, without destroying the coherence properties of the latter. 

We start by defining three energy scales, viz., a = naT, (J"^ + A^)^/^ and A, where A is defined as 



A^ = 



A 



*2 



/27rT 
r(l - 2a) \ A*~ 



2a 



(5.34) 



Of these 3 energy scales, a is the smallest. It is now convenient to define (by analogy with the biased spin- boson 
modelEl): 



7fl,s(a/3,>77 A^) = 2ap 



A2 



J^ + Aj 



jsiajs, J, Aj) = 2ap - - jr.. 



af3 



J^ + Aj 



(2^2 + ^2 ) 



^J^ + A| - 



,2 a4 
J 



ajAj 



(5.35) 



j2 



with a 13 



■KOtpT and Ap defined as in Eq. 
2 „, 2a« and - j'^. 



2apAl/J 



Is 



(5.34) with A replaced by A^. In the limits considered in this section, 



1. Equivalent Spin Case 



For this case, defined by T ^ and A/a ^ T for both spins, the equation for the roots of the denominator in 
^TiT2(A) decouples into 2 equations of polynomial degree 3, which we can handle by perturbation in a = ■naT. The 
resulting forms for the "FM- correlated matrix elements" PT,'[{t) and Pii{t) are then found to be (cf. App. E): 
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-Fat 



Cae"^ * cos(^<) + C4e"^ ' sin(/it) 



(5.36) 



Pxi(t) = ^--C2e 



+ C3e ^*cos(^i) + C4e^^'sin(^<) + 
Cg-e--^'* cos(j^,t) + ^7-6-''* siii(i/,t) 



(5.37) 



where Tr = 777,^(0, J, 2A2), T = ■^,{a,J,2A?), = vl{a. J,2A^), and A± is_defined by Eq. ( |5.25| ). If we now 



introduce J'l^ = ^JA^ and Aj_ = 4A"^A±, then we have ^r± = 7i?,s(a, J'lp, Aj_), 7± = 7s(a, J'lp, Aj.) and the 



frequencies k 



(a, i/ip, Aj_). The exphcit expressions for the 3 characteristic frequencies of the system are then: 



J2 + 



(J72 + 2A2)2 



j2- 



(5.38) 



a^A^ 



.2\2 



AAa 



16A: 



'A2 



(5.39) 



The frequency fi pertains to the combined system, whereas the fre quen cy pair iy± desc ribes correlated oscillations of 
the individual spins ("beats"). The various constant C^'s in Eq. ( ^.3(j ) and Eq. ( 5.37 ) are complicated functions of 
the frequencies and decay rates and given in Appendix E. 

On the other hand, in this equivalent spin limit, the ATM correlated probabilities Puit) and Pi^it) have a much 
simpler form: 



Pn{t) = Pi^{t)=A+ + C2e 



€36 '"*cos(^i) - 6*46 '"*sin(^i) 



(5.40) 



Let us now try to understand some of the details of this result. The appearance of oscillations in the density matrix 
at 3 different frequencies reflects the fact that the 2 spins are neither in a completely correlated nor in a completely 
uncorrelated state. The total dynamics then results from a compromise between global and individual dynamics. The 
global dynamics can be eliminated if one considers the probability of a spin being in a given state, independently of 
the state of the other. For example, the quantity [-Pfi(i) + ^'Ti(^)] describes the dynamics of the first spin, with the 
second one "avera ged o ut" . Notice, incidentally, that it is not possible to recover the oscillation frequencies of our 
"toy-model", Eq. ( |5.38| ) and ( ^.39D , in the hmit a = since paths with overlapping blips do not contribute to Ptit2 
in the uncorrelated phase. The lower limit on a is set by ~ A/a. 

To leading order in T/J and A/ J, the decay rates Tji and 7^,^ are given by 77?, i = Tii/2 ~ 2TTaTA'^/J', which 
simply corresponds to the hmit aT <C JT" (always with ^7 A) of the individual decay rate given in Eq. ( ^.30 ). 
The same is true of the decay rates F and 7^, the decay rates associated to the oscillatory terms. However, they are 
strongly corrected by a term of order 0{{J / A)"^), so that the oscillations are damped out very quickly compared to 
the general decay of the probability function. 

We illustrate this case with Fig. (|To|), showing a special case of the Fourier transform of P|f(t), namely 
/m[— AP|f (iA)], Without the environment, the oscillation frequencies would show up as divergences in this function. 
With the environment, the divergences are smeared out by the decay rates F and 7±. At this scale, the differences 
between /i and v± are clear, but not the difference between and v^. Taking a larger value of a (resulting in a 
higher a — naT) increases the decay rates associated with the oscillations and the width of the individual resonances 
becomes large enough for the two to overlap. 



2. Overdamped plus Underdamped Case 

We can also imagine two different systems labelled by 1 and 2, such that A2/a2 > T whereas Ai/ai < T. In a 
single spin analysis, the spin-1 would be underdamped with the spin-2 overdamped. The limit 3> Ai, A2 assuring 
the non-overlap of the blips is implicit. 

To establish the dynamical behaviour, let us first look at the two denominators involved in the Laplace transform. 
In Appendix F, we show that the global denominator (A-|- gi + g2) yields decay rates and oscillation frequencies given 



by F2 -I- "fR, F2 -I- 7 and /i , where F2 is the decay rate of the overdamped spin, Eq. (5.3C) and 7^,, 7 and /i 



expressed with respect to Eq. (5.35) as 
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^ = -f,{a,J,K\) (5.41) 

where a = ai — T2/2. 

To treat the second denominator [(A + gi)(A + 52) — we simply drop the terms /ii (A) /12(A). We can do 

this because hi ~ J /T, whether the spin is over- or underdamped. Thus as far as the decay rates and oscillation 
frequencies are concerned, these terms will only be a small correction of order (J'/2T)^(Ai/J')^ ^ 1. We are then 
left with the decay rates and oscillation frequencies of a single over- or underdamped spin in a static bias . These 



are r2, as give n by Eq. ( |5.30| ), 71, ^ = 7^,^5(01, J', Ai), 71 = 7^(01, j7,Ai) and the frequency z/f = i^^(ai, J', Ai), 
defined in Eq. ( |3.35 ). Notice however, that the correlations brought by hi, /12 in the numerators of PriTii'^) must be 
kept to determine the weight of the relaxation and oscillation terms as well as the long time limit of the probabilities. 

With the decay rates and oscillation frequency identified, it is straightforward to perform the inverse Laplace 
transform to obtain the probability of occupation P-\^{t) as 

Pn(t) = A_ + C2e-^^-^« + Cse-^^-^ cos(/2i) + de-^^'^ sin(/2t) 

+ Cge-^^* + Cee^'^^-^* -f Cye"''!* cos(z^it) + Cge"'^!* sin{vit) (5.42) 

Again, oscillations associated to the global and the individual dynamics of the spins are present in Notice 
that the corrections to the frequency ji are "tunable" if one has very good control over all the parameters. In fact, 
with, r2 = fli, the frequency ji^ is simply equal to + A^. On the other hand, once again most of the effect of the 
environment is felt as a damping of the oscillatory terms, much faster than the the global decay (with the oscillations 
of frequency /2 being damped faster than those of frequency vi). 

This case is illustrated of Fig. (pT|). The coherence peaks are still present, but the line- width caused by the 
environment impedes their resolution. 



VI. CONCLUSION 



We have spoken very little in this article about the possible physical applications of our results, preferring first to 
gain a general understanding of how the model behaves. We have in fact already used these results to do a number of 
detailed calculations for 2 particular problems. The first concerns the dynamics of a pair of nanomagnets, each coupled 
to a background conducting or semiconducting substrate. To properly treat this we need to add to what we have done 
here the presence of an external bias acting on each nanomagnet. It is straightforward to make this generalisation, 
and also to calculate from microscopic theory what are the parameters entering the PISCES Hamiltonian as functions 
of nanomagnet size, distance between nanomagnets, electron density, etc.; and one may also incorporate the effect 
of nuclear spins and dipolar interactions. This problem is not only of academic interest; it is also relevant to future 
information storage technology. 

The second problem we have looked at describes a fairly general class of physical systems in which we have a 
macroscopic measuring apparatus coupled to a central "measured" system; and in which at low energies the dynamics 
of both truncates to 2-level systems. The purpose of this work is to go beyond the usual idealised models of the 
measurement process, and see how the 3 important partners in the measurement (system, apparatus, and environment) 
work together. Again one finds that the PISCES model (with added biases) can be used to derive the full dynamics, 
with some rather surprising results. Both of these investigations will be the subject of detailed future papers; and 
they confirm that the results we have here can be used directly to give a realistic description of the full quantum 
dynamics of coupled macroscopic systems (as well as a large number of microscopic ones) . 
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APPENDIX A: REDUCED DENSITY MATRIX OF THE TWO-SPIN SYSTEM 



In this appendix, we derive the form of the reduced density matrix of the combined two-spin system. We will 
assume the system starts in the polarised state | tT)7 so both 7710 and 7720 must be +1. Now consider the density 
matrix element Ptit2 ,T\ (t) i i-G-, the probability -Prir2(0 that at time t the system is in the state T1T2 where Ti and 
T2 represent either t or |. Using the decomposition of the influence functional into individual functionals and the 
interaction F12, we express the complete summation for the probabilities as 

Fi^^ ({i. }, {Ci,}, {^1,})^^^^ {{U2k}, {C2k}. {mk})Gn,n. {{h,} , {u^u} . {Ci,}, {'71, }, {Csfc}, {mk}) (Al) 

We first explain each term of this expression. The factors d)^p are different for each probability calculated and reflect 
the constraint on the last sojourn of each paths. If we are interested in calculating (f), then at time t, the two paths 
must be in the state []]] after having performed 2ni and 2n2 jumps. We thus have the constraints 77i„j = 772712 = +1 
and dji^ = cfy^ = 0. However, for Pii{t), the paths must end in the state [H]- This means that at least one blip must 
be present in the paths considered. They will thus be composed of 721 + 1 and 7x2 + 1 blips, both 711 and n2 ranging 
from zero to infinity, with the constraints 77irii+i = 112112+1 = ^1 and dj^' = = 2. Finally, of course, the paths for 
P'\l{t) and Pi-\{t) will be formed by a mix of these two cases. For P|j, the constraints are 771,^ = +1 and 772712+1 = ^1 
with (f^^ = and d^^^ = 2, while for the constraints are ?7irii+i = —1 and 772ri2 — +1 a-nd d^^^ ~ 2, ~ 0, 
The time integrations are defined as 

D{t2n} - / dt2n / dt2n-l ' ' ' / ^^2 / dti (A2) 

Jo Jo Jo Jo 

ie., an integration over the times at which the jumps occur. Notice that the blips of a given path are well ordered, 
ie., they occur successively; however, nothing specifies the order of the blips of one path with respect to those of the 
other, and so all the different different configurations must be included. For example, consider a path relevant to the 
calculation of P]]{t) with ni = n2 = 1; there are four transition times ti, t2, ui and U2, which can be arranged in 6 
different configurations, as represented in Fig. ^2|, always with the restriction ti < t2 < t and t^i < 7^2 < t- In general, 
there will be 

2r7i + 2772 \ (2^1 + 2772)! . 

2771 )~ 277l!2772! ^ ' 

such configurations. This is simply a reflection of the fact that we cannot describe the two paths with respect to a 
single time where the all the jumps would be well ordered 

The terms Fr!^ {{tj\^ {Caj}, {Vaj}) are the single-spin functionals which can be decomposed as 

where 



pia;se) 



exp 



2 

■— Q2°'\t2j - t2j-i) 

TT ^ — ' 

i=i 



(A5) 



P(a:h-b) _ 



exp 



4 



j>k 



(A6) 



exp 



2 

i^VlT -77 tX^"^ 
TT ^ — ' 

j>k 



(A7) 
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in which the kernel matrices are 



^ik = Q2'\^2j — t2k-l) + Q2^\t2j-l — t2k) — (^2j — t2k) — Q2"''(^2j"-1 — t2k-l) 



(A8) 



^jfc'' ~ Ql ~ t2k+l) + Qi '{t2j-l ~ t2k) — Q\ '{t2j — t2k) ~ Qi ' [^2]-! ~ i2fe+l) 

and are defined in terms of the correlation functions 



duj 'I^^ [ 1 _ coswt ] coth(c<j/2T) . 



(A9) 



(AlO) 



CLLU T\ — sm uji 



(All) 



The correlation Q^"\t) refers to the reactive (phase) term ^^(i) in ( 4.22 ), and Q'^'{t) refers to the dissipative term 

As discussed at great length by Leggett et al., ( |Al[ ) can be analysed in many regimes using the "dilute blip" 

approximation, in which the interactions A^-^'' are dropped entirely and only the term f'i"'*'^'* and the term in A^j^'' 
corresponding to the interaction of a blip with the sojourn immediately preceding it are kept. Thus the single-spin 
influence functionals -FtI"'* in ( Al ) can be evaluated without worrying about the precise ordering of the charges in the 
configuration. 



Now let us turn to the interaction hmctional F12 in (4.23). One finds that F12 is a sum over all possible double 
paths integrals of the form 

[ D{t2nA f D{u2n,} J] ^ G„,„, ({* i, } , } , {Clj } , {^7l, }, {C2j }, fej }) (Al2) 

{Cl3,»?lj} {C23,»)2j} 

with weightings (— iAi/2)^"i (— iA2/2)^"^ ( here we have chosen a set of path corresponding to the calculation of 



P||(t); (c.f., Eq (Al)). The integrand G„j„2 has the form 

G„,„,(R) = Gi''-„t)(R)G(''-«,)(R) 

where 



(A13) 



G(''^-;,^,)(R) = exp 



<Zoi9o2 



TT 



ci]C2k^jk{'R') 



(A14) 



in which the blip-blip interaction matrix is 
l(l2)/■D^ _ r^(i2)^T, 1. „, h , ^(12) 



A)k>m = [Ql'^'iR, \t2j ~ U2k-l\) + 0^"'(R, \t2j-l ~ U2k\) - Q^2^\'R, \t2,-l - U2k-l\) ~ Q^"' (R, \t2j - U2k 



(A15) 



Notice the difference from (A8). The Gnin2 C^) term is 



G(f-«^)(R) = exp 



.'Z01902 



TT 



rii n2 



^^Ajf (R)(Ci,?72/c+r7i,C2fe) 



exp 



j k 
rii 



-*i^z.(R) 



Clj?72p(^2j - t2j-l) H- ^ C,2k'nip{u2k - U2k-l) 
j=l k=l 



(A16) 



where 



Ajf )(R) = QY^'CR, \t2j - U2k-l\) + Qi"'(R, \t2j~l - U2k\) - 1^2,-1 - U2k-l\) - Ql"'(R, \t2j - U2k\) 



1(12), 



l(12) 



(12), 



1(12), 



(A17) 
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Kzz is the direct interaction and couples blips and sojourns belonging to different systems, but only while they overlap 
in time. A term directly coupling blips with earlier sojourns would come from retardation effects, and these are ignored 
here. This restriction is indicated by the index P appearing in the charge of the sojourn. The correlation functions 
have the obvious definitions 

Q(i^Ht) = r dcj sin^t (A18) 

Jo ^ 



Q^^^\t)= f dio -^1^^^'^) [1 _ cosu;t]coth(^/2r) (A19) 
Jo 

We notice that there is no term in ( [A13| ) corresponding to the term F!i^^ in (A4), i.e., no "self-energy" term in the 
interaction functional corresponding to the the blip self-energy in the single-spin functional Fj] G„j„2 contains only 
spin-spin interaction terms. This point is crucial - if we try to make an approximation analogous to the non- interacting 
blip approximation, so useful for the single spin-boson problem, we would be throwing away all spin-spin interaction 
terms! Thus the PISCES problem is intrinsically more difficult than the single-spin problem, since we must in some 
way deal with these interactions. 

Consider now the structure of the interaction influence functional. We need only consider three explicit 
configurations- the main characteristics of the analysis depend on whether or not the blips of a path overlap the 
sojourns of the other (well separated blips and sojourns give standard results). First, if the blip and sojourn are not 
overlapping, the interaction is identical to what was found before. An example of this case is the interaction of the 



first spin-2 blip (charge C21) with the first spin-1 sojourn (charge t^iq) in Fig. 12-a The result is: 



exp z 7?io C21 [ (R, - io) + Qi'"^ (R, u2-h)- g^^'^ (R, u2-h)- Q^''^ (R, ui - ii ) ] (A20) 

However, if a blip occurs entirely within the time interval of a sojourn, complications arise. As an example, consider 
the interaction of the first spin-1 blip (charge Cii) with the first spin-2 sojourn (charge 7720) in Fig- P^b The time 
integration is: 



r 11 

/ dT I ds sitllo{t — s) = —[t2 — ti) -[s\nuj{t2 — to) — sm-uj{ti — to)] 

Jti Jto ^ ^ 



(A21) 



which gives the result 

exp t Cii V20 [ ( (R, h - to) + Q[^^^ (R, t2 -to))- {Kzz (R) + e(R))(i2 - ti) ] • (A22) 

where e(R) is given by 

,(R)^» Trf^M^. (A23) 



This function is discussed in the main text. 

The 3'"'' case involves overlapping blips, and the total interaction appears as well. If the blips are overlapped as in 
Fig. 12 -b it is convenient to consider together the Cn — 7720 and the C21 — ?7ii interactions. The corresponding factor 
in the influence functional is then 

exp<ii772o[Ql^^^(r,i2 - ui) + (3j^^^(R,ti - to) - (3i^^^(R,t2 - to)] x exp<2i??iiQl^^'(R,U2 - ^2) 
X expiJ^[Cii772o(ui - ti) + C2i??ii(u2 - ^2)] (A24) 

Finally, if a blip completely overlaps another blip, as in Fig. |l2[-c, the result of the interaction with the two neigh- 
bouring sojourns is 

exp-|-iJ'[Cii772o(wi - ti) + Cll?72l(t2 - U2)] 
expzCiir72o[Ql''HR,ti - to) + Q1''^(R,<2 - «i) - Ol'''(R,<2 - to)] - iCiimiQi^\'^M - U2)] (A25) 

All other cases in Fig. |l2|are similar to these, provided one interchanges the times tj and Uk- 
This summarises the construction of the influence functional for any conflguration. 
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APPENDIX B: POOR MAN'S SCALING FOR THE OHMIC PISCES PROBLEM 



We start from the reduced partition function of the Ohmic problem, expressed as: 



Z/Zenv = / D[qi]D[q2\e 



-So[gi]-So[<j2]+J(R) pj"^ dsqi(s)q2(s) 



I fi/T rl/T 

exp — dsj^ ds'[T]iqi{s)qiis')+rj2q2{s)q2is')+Vl2{'R)iqi{^^^^ (Bl) 

where Zenv is the partition function of the environment, 5*0 [qa] is the action of the a-th two-level system, and the 
remaining terms represent their interaction with the oscillator bath. Notice that the total adiabatic interaction 
J^(R)qi{s)q2{s) is used, so that all interactions resulting from static configurations are explicitly separated from those 
coming from changes in the configurations (ie., transitions). Each term in the equation can easily be interpreted with 
respect to the single spin-boson problem. 

The partition function is evaluated using the usual instanton method. The trajectories are defined by a series of 
tunneling jumps (the instantons), of duration smaller than Tc, taking place at times tj and connecting the different 
configurations of the system. We refer to the different configurations by the label zi, I = 1,2,3,4 such that z = 
{{11}: {11}: {it})- These are connected by the tunneling matrix elements A(fc, I) — A{1, k); here A(l, 1) = Ai 
corresponds to transitions between zi or Z2 with z^ or Z4, A(2,2) — A2 connects 21,24 to 22,23 and we introduce 
A(l,3) = Ap and A(2,4) — Aaf connecting the ferro- and antiferromagnetic configurations 21 — 23 and 22 — 24 
respectively. Even if no such transitions are present in the original problem, they will be generated through the 
renormalisation procedure. We also define the dimensionless fugacities y{k,l) = TcA{k,l). To each configuration 
is associated a static energy e;, and a corresponding dimensionless energy Mi = T^ei- For the PISCES problem 
considered here, Mi — M^ ~ TcJ /2 and M2 = Mi = —TcJ/^. 

The boundary conditions on the trajectories are such that a path starting in a configuration 2; at t must return 
to this configuration at the time t — 1/T. Paths starting from the 4 different configurations must be considered. The 
time integrals in Eq.(Bl) can now be done, and the partition function is 



Z/Z, 



n=0 



n 2i 1=1 



^'^ dt, 

y(M-l)/ —QiU 

'''c 



Tc) 



exp 



E 

i=a 



ti - ti^ 



Mi + U{{ti}) 



(B2) 



Just as for the dynamical behaviour of the PISCES problem, nothing specifies the ordering of the instantons and we 
must sum over all the possible set of configuration {2/}. The time integral accounts for all the possible positions of 
the instantons, with the O- function insuring that none are within Tc of each other. The effect of the bath is contained 
in the interaction matrix defined as 



U{{ti})^Y.^n 

l>3 



U — U-i 



(B3) 



with Ki j representing the interaction of the I instanton with all those at earlier times. In terms of the parameters 
of the PISCES problem, these are Kfj^2 = a/32, K/s^p = a/3 + 0:12, Kp,AF = otp — ai2, and, Kp^p = ai + a2 + 2ai2, 
Kaf,af = ai+a2 — 2ai2, Kp^AF = ai — a2- This interaction is logarithmic, which is what allows the use of the poor 
man's renormalisation group. If one suppresses the energy difference between configurations, this partition function 
is completely equivalent to the one treated by CardycJ. In the present case where symmetry between the levels is 
broken, the poor man's renormalisation procedure requires the inequalities y{k,l) <C 1, T ^ fig, and <^ Qq. If 
we also assume the 2 spins are identical, ie., that Aj-^ Ai — A2 and a = ai = 0,2-, then the partition function is 
equivalent to that derived by Chakravarty and Hirschllll for the 2-impurity Anderson model, and we may simply take 
over their scaling equations, quoted in the text. 



APPENDIX C: DYNAMICS IN THE LOCKED PHASE 



When the indirect coupling is the largest energy scale, the two spins will tend to tunnel simultaneously. The 
time interval spent in a "unlocked" state ( {Ti},{iT} if the coupling is ferromagnetic (or {tT},{ii} if the couphng 
is antiferromagnetic ) will be ~ In terms of our path integral for the density matrix, this means that the 

dominant paths will be those where the blips are overlapping. This minimises the overlap between blips and sojourns, 
thus cancelling the fast varying factor expjj7(tj — u^) in the influence functional which would otherwise give this 
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configuration a vanishing contribution to the path integral. Thus, we can restrict the summation to be only over the 
paths where the beginning (end) of one blip is within J"^ of the beginning (end) of the second blip. Furthermore, the 
paths <7i(t) and qi(T') are nearly identical to the paths q2(j') and 52 (t') respectively, and this brings a restriction on 
the values that the charges of the blips and sojourn can take. In the ferromagnetic case, they must be equal. That is 
Cij = = il and r\\j = r\2j = il- However, in the antiferromagnetic case, they must be opposite (C^\j = —C2j = il 
and r]ij = -772^ = ±)- 



With these paths, the time integration over the end-points of the blips, Eq. (Al) is greatly simplified. For strong 
bias, such that is the shortest time scale of the problem, we can set tj — Uk = tj — tk provided that j ^ k. In the 

f I2I fl2) (12") (12") 

influence functional, we can also set Qr'ihj - U2j) = Q\ i^) = and Qr'{t2j - U2j) = Q\ ' {<d) = since they 
will be of order ll\J\. We clearly cannot do the same with the factor expij'{tk ~ Uk) but it is now easy to integrate 
over it. Always in the limit of strong coupling, the integration over the the endpoints of the n*'' pair of overlapping 
blips with charge placed between sojourns of charge rjn-i and rjn will be 

r^/J r^iJ , / +1 \2 1 

4 lim / dwe^*""-!'^-^" / du'e^*""'^-^" -r-r T^— (CI) 

^/J^^Jo Jo \AJ\J Vn-lCiVn 

u and u' being short times over which the blips are not overlapping. The signs ± correspond to ferro- or antiferro- 
magnetic coupling respectively and the factor 4 comes from the four possible configuration of two overlapping blips. 



In Eq. (CI), we have neglected a factor of 0(1) coming from the terms e*''"''" — 1). Taking into account the {n -\- 1)*'* 
blip brings a equivalent factor, but with charges f7r!Cn+i^»i+2 so that for a chain of n overlapping blips, the total 
contribution from the time integration of the overlap is 

' ' (C2) 



This shows that we now have a new effective tunneling matrix element equal to Ac = AiA2/|i/| (compare Fig. (|^); 
there, the problem is set up so th at Jp = i//2). Furthermore, within this approximation, the interaction of two 



overlapping blips, as given by Eq. (|A15|) becomes 

A,fe(R) -> 2Qi'^\t2k - t2k-i) (C3) 

which is simply a contribution to the self-energy of each blips. Similarly, the interaction of the blip of one spin with 
the sojourn of the other spin becomes also simply adds to the contribution of the individual spins. 

Therefore, the whole system behaves like a single spin-boson system to which is associated a tunneling matrix 
element Ac and coupled to the environment through a new coefficient ac — ai + a2 + 2ai2 for ferromagnetic coupling 
(since CijC^j = +1) and ac = ai + ^2 — 2q:i2 for antiferromagnetic (in this case CijC2j = ^1)- 

Now the dissipative dynamics for this system is very well known, and all known results can be taken over directly. 
Assuming a ferromagnetic coupling, so that the system oscillates between the state {tt} and {H}, it is straightforward 



to use the dilute-blip approximation and to obtain Eq. (5.6), the Laplace transform of P||(t). The expression for 



P|l(A) will be equivalent, but with the -I- sign replaced by a — sign (coming from the constraint 770 = 0, 77„ = —1. 



with /(A) given by Eq. (5.1C) 



Notice that if we had an antiferromagnetically coupled system, starting at i = in the position {ti} or {jt}i the 
structure of the summation is completely similar, one merely needs to replace P|| and Pn by P^i and P[|. The case 
of antiferromagnetic coupling, but starting configuration {tt} is discussed in the main body of the text. 



APPENDIX D: DYNAMICS IN THE CORRELATED RELAXATION PHASE 



As long as T 3> i/, the renormalisation group analysis shows that there is no great tendency of the 2 spins to lock 
together. In terms of the path integral, this means that the oscillating factors exp{i(^ajill3k>Jitj — Uk)) do not reduce 
the weight of paths with overlapping blips and sojourns as efficiently as in the locked phase. On the other hand, in 
the limit J ^ A^, if the blip of one path overlaps with a blip of the other, a factor A1A2/J' <C A^ still appears in 
the path integral, which then tends to reduce drastically such overlaps. If furthermore, the blips on the individual 
paths are naturally dilute, then the occurrence of overlapping blips will be extremely rare and we can neglect them 



25 



completely. The paths summation can then be done by using only paths where a blip of one path only overlaps in 
time with a sojourn of the other path. The two essential conditions of this approximation are thus T 3> JT", to insure 
that we are not in the locked phase, and J ^ Ajs so as to neglect paths with overlapping blips. In addition the blips 
of each individual spin must be in a dilute statead. If a/j > 1, this is so for all values of J and T. On the other hand, 
if ^ 1, then the blips are dilute provided that 



dX 



A *2 



A=0 



i2n apTf + 



(Dl) 



Notice that as long as ^ A^, the ratio apT j can be arbitrary, thus allowing the occurrence of the mutual 
coherence phase. 

We can now proceed to explain in detail the path summations that lead to the expression given in section V-C. 
Each paths are formed of n\ blips of the first path, and rii blips of the second path, such that no blips are overlapping. 
This approximation corresponds to the dilute-blip approximation for the single-spin influence functional and means 
that we keep only the direct interaction (J7 — Kzz + e) in the interspin functional G„j„2. All the other terms vanish 
by the usual arguments of the dilute-blip approximation. 

Blips of a given type are well ordered amongst themselves. What changes from one configuration to the other 



is the ordering between blips of different type. Therefore, the summation over the {Cij} and {C2/c} in Eq- (Al) is 
independent of the configuration and can be performed immediately. It brings a factor 

rii 

2"! W COS [mi-iQiHh] - t2j-i) - mpJihj - t2j-i) 

2"=^]^ cos ri2k^iQf\u2k ~ U2k~i) - VipJ{'>J-2k - U2k-i) (D2) 



3 = i 

X 

fc=l 



where rjip and ri2P are the charges of the sojourns overlapping the blips ^2*; and respectively. 

With the charges of the blips removed from the problem, we will refer to the blips by the charge of the sojourns 
immediately preceding them. The fact that the blips are not overlapeing then allows the use of the Laplace transform 
to perform the summation, just as in the single spin-boson problemld. To keep track of the different indices involved 
in the summation, it is convenient to define a 2 x 2 matrix g(")(A) with components (g'-"')^ as|^ 



(D3) 



with the definition 

.2 



A 1'°° 

(.g(")(A))^ ---f / dte-'''e-'^- ^'^cos[Tjt~aQ["\t)] (D4) 
^A Jo 

The function (g^"^)^ represents the contribution of a blip of type-a. The index a refers to the charge of the sojourn 
on the path a immediately preceding the blip while the index r refers to the charge of the sojourn overlapping the 
blip. 

The contribution of the blips to a given path will now be expressed as a product of the components of g'^")(A). What 
remains to be done is the sum over the configurations and the set of the and {?72fc}- This must be performed 

simultaneously since the summation over the rj gives a different result from one configuration to the other. 

We will give a detailed explanation of the procedure to calculate P^^{t) and simply state what needs to be done 
differently in calculating the three other probabilities. The basic building blocks of the summation are the chains of 
coupled "clusters" of blips. The clusters are defined as the ensemble of successive blips of the same type a, successive 
meaning that they are all overlapping the same sojourn of charge ?7;3,_p. This corresponds to a single spin-boson 
system evolving in a static bias rjp.pj. All the different clusters are then linked together by the charge of the sojourn 
preceding the first blip of a given cluster. P]]{\) then consists in a summation over all the possible chains of clusters. 



^We thank E. Mueller for pointing out this matrix method to us. The complete path summation was first done using a more 
complicated combinatorial method. 
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each cluster containing all the possible number of successive blips. This is represented schematically on Fig. (|l|) and 
Fig. dll). 

The complete summation over a cluster is easy since it is equivalent to the spin boson problem. The total con- 
tribution of a cluster of type a, whose preceding sojourn has a charge a and overlapping a sojourn of charge r is 
then 



{a,} 



1 



l-(5("))^-(.9M)r 



(D5) 



Notice that the expression (^("^)+ — (^("))~ is nothing but the function P(A) = P^iX) — PiW obtained for the single 
biased spin-boson system in a static bias rjT" within the dilute-blip approximation. 

The summation of the chains of clusters is now straightforward. Let us refer to the contribution of a chain of 
clusters beginning with a cluster of type a and ending with a cluster of type (3, with n clusters of type a as C^'^ . The 
summation over the charges of the sojourn is included in this notation. The complete probability summation can be 
expressed as 



Pi 



TT 



1 1 

A + A 



OO -| OO OO -| oo 



71=1 



n=l 



(D6) 



n=l 



n=l 



where the factor 1/A comes from the definition of the Laplace transformation. Let us consider in detail a particular 
contribution P"~^ = ^ C"^"- It is composed of a sum over a product of 2n — 2 successive clusters (n a-clusters and 
n — 2 /3-clusters). No clusters of the same type are adjacent. Before the summation over the charges, the product 
is thus of the form (g''"'')^^' ...(5^"')tJ (5^'^)tJ+i ■■■(5^"'')T|r-2 ■ However, since they are successive, there is the 
restriction crj+i = tj. It is then possible to write 



(D7) 



n=0 



3 = ^^,} 



where for each n, there is a restriction (72n+i — 1- This is then nothing but a particular element of the multiplication 
of two matrices: 



p(a-/3) ^ (^("))+ g(/3)(l _ g(")g(/3))-l 

The similar expression for a chain where the two limiting clusters are of the same type is 



(D8) 



n=lj=l {aj} 



The complete probability P||(A) is thus given by 



^Tt(A) = \ + \i9^'Y+ [(1 - g^^^g^^^)-^] I + \i~9^'Y+ [(1 - i^'^e^y 

(^(1))+ [g(2)(l-g(l)g(2))-ll+ + Cg(2))+ [gW(l-g(2)gW)-l 



(D9) 



(DIO) 



The other probabilitie s are then obtained in a straightforward manner, they simply correspond to taking different 
matrix elements in Eq. (DIO). We obtain: 



5-(i)g(2)(i _ g(i)a.(2)^-i 



(Dll) 
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Notice that the matrices now involve [• • -jl , since the final sojourns must have a charge of —1. There is also the factor 
1/A missing, since there must be at least one blip in each path, to go from the state | tt) to | ||). 
Finally, for the last two remaining probabilities. 



+ A' 



5-(l),3-(2)^-l 



(l_gU)g 



g(l)(l_g(2)g(l))-l 



(D12) 



with Pxi(A) simply related to P^iX) by the substitution (1) <^ (2). 

These are obviously fairly complicated expressions, despite their compactness. However, the matrices g^"-' are 2x2 
so that it is relatively easy to expand. Notice first that the denominator of Ptit2Wi which determines to poles in A 
and consequently the dynamics is given as 



1 -rr(g(i'g(2)) +i:»et(g(i)g(2)) 



(D13) 



This denominator reduces to a product of polynomials of degree 2 in A and it is then possible to simplify PT1T2W 
further. 

At this point, we can revert to the n otation of Leggett et al. and use the functions g/3(A) and hfj{X defined by Eq. 
( Ft5|) to simplify P^^iX) to Eq. (^TgI) of section V. 
Similarly, we obtain for Pii(A) 



91 1 



91 



92 1 



h2 
92 



1 



1 



A X + gi+ 92 



51 ( 1 + — ) + -92 fl + — 
91 ) \ 92 



1 



A2 + A(gi +.92) +31.92(1-^) 



-.9192 



1 



hxh2 



9192 



1 



A A2 + A(5i +. 92)+ 5152(1 -ff) 



While the last two remaining probabilities simplify to 
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~ 4 



<7i 1 + - 









_h2 




+ 92 
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1 



1 



A A + .91 + .92 



1 



A2 + X{gi +92) +9192(1 



hih2 ' 
31 32 ' 



(D14) 



(D15) 



At (A) = 7 



51(1 + ^ 

91 

9i{l + ^ 

91 



h2 

52 1 + — 

92 

52(1 + ^ 

52 



1 



1 



A A + 51 + 52 



A2 + A(5i +52) +.9152(1 -ft) 



(D16) 



For identical systems, P^i = Pn, as should be expected, but the form of these probabilities also become quite simple, 
setting 5 = 51 = 52 and h — hi = h2, the probabilities are 



(Identical systems) 



(D17) 



APPENDIX E: MUTUAL COHERENCE REGIME: IDENTICAL SYSTEMS 



In this appendix, we discuss the more tedious details pertaining to the mutual coherence case. We start with 
the denominators in Pr-^r2W^ whi ch yi elds th e osci llation freq uenc ies and decay rates and then give the form of the 
coefficients Cij appearing in Eq. ( 5.36 ), Eq. ( 5.37 ) and Eq. ( ^.4C| ). All the following expression are of course quite 
messy, although straightforward to obtain. Furthermore, they are obtained by working to lowest order in a = naT. 



28 



It is nevertheless important to get an idea of the value of relaxation rates and oscillation frequencies as well as the 

relative weight of the different terms in the overall probability function. 

For two equivalent spins, the denominators of the Laplace transforms Prir2(A), Eq. ( ^.17 ) and Eq. (5.18) are now 



A + 25(A) (El) 

and 

i\ + giX) + hiX))iX + giX)-h{\) (E2) 
With the appropriate forms of g{X) and h{X) inserted, the equations for the poles of the Laplace transform are then 

X^ + AaX^ + XiAa^ + + A^) + 2aA^ (E3) 

+ AaX^ + A(4a2 + + A^) + A^{2a ± J) (E4) 

The structure of these equations is completely similar to the equation obtained by Weiss and Wollensack for the single 
biased spin-boson systemtll. The equations have two imaginary roots if T > A/a, which corresponds to oscillating 
terms in PT--^r2{t)- The form of the decay rates and oscillation frequency can b e exp ressed easily in terms of the 



relaxation rates and oscillation frequency of the single spin-boson problem, Eq. ( 5.35 ) and yield the values cited in 
the text. 

We can now discuss the values of the coefficients Cy . The coefficients C2, C3 and C4 are related to the dynamics 
of the two systems together. C2 is the coefficient of the pure relaxation term, with value 

AV2 / 8aA+\ 

= (Twtt^ V " ^ J ^""'^ 

In the limit J' ^ A this coefficient is of order 0(1), with limiting value 

C2 ~ A_ iJ » A) (E6) 

which makes the connection with the overdamped relaxation of the correlated phase. The coefficients C3 and C4 
corresponds to oscillating terms with frequency /i. It simpler to write them directly as a sum of cosine and sine since 
it avoids the complication of introducing a phase in the oscillations. We obtain 

^ A^ r^ + A^^ + 8aA+(rfi-2r) 

The main point about these coefficients is their smallness. In the limit J ^ A^ they are 

C, ~ ^ (E9) 

(El„, 

We now get to the coefficients pertaining to the individual dynamics of the two spins. In this case, the Laplace 
transform of i^rir2(0 includes terms with denominator of the form (A -I- ^(A))^ — /i^(A), giving rise to very complex 
coefficients. There is the apparition of two individual relaxation rate 7_r±, with a coefficient CJ^, where the index 
T = ±1 refers to -Pfi(t) (r = +1) and Pii(t) (r = — 1). Its expression is 



In the hmit of interest to us, J ^ A, Cl^ is given by 
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Which is of order 0(1), as is expected since we are close to the totally overdamped limit. 

The remaining terms are the oscillations of frequency i'± and damped at rates 7±, with coefficients Cg^ and Cf^, 
where r refers once again to Pf| or P|j. Of course the exact form of these coefficients is fairly complex: 

^^''^ ^ 2^^^ ' + ^:^2^ - N^D^2)) (E13) 

^^'^ ^ ^^J^ + ^^'^^ + "^|T^(^3^^? + ^2 )) (E14) 

The different terms appearing in the numerators and denominators are given by, 

= [--/a + 8aA+][{2a - 7^)2 - + J^] - 2vl{2a - 7,) 

~ -SaKlA+ (E15) 
N!^ = 2v„{2a - 7„)h7a + 8aA+] + i/,[(2a - 7,)^ - ^ 

^ -JKI (EI6) 

iVa" = -7.[(2a - 7,)^ - i.^ _ 4^2 tanh2(^/2r)] + 2!y2(2a - 7,) 

- ~laJ^ - 8J^aA+ (E17) 

Nl - -i.,[(2a - 7^)2 - j/2 _ V tanh'(j7/2T)] - 27,i^,(2a - 7,) 

- (E18) 

The dominant term in then N^fj^ and the difference between Af4+ and N4^_ can be neglected at this level. 
The remaining terms are given by 

Dl = [(7- - 1+f + ^-^ - ^l][{lR+ - l<y){lR~ - la) - vl\ - '^^lil-a - la)[lR+ + IR- " 27,] 

^ -AK'j-^{A_^~ A,) (E19) 

D2 = '2Va{l~a ~ lcr)[iriR+ " 7<T)(7i?- " la) ~ + V„[-iR+ + 7ij_ - 27,,] [(7^ - 7+)^ + v'i^ - vl] 

-Sajh'^{A-a- Aa) (E20) 

\D''? = [(7fl.+ - 7.)' + ^IWilR- - 7.)' + ^IWl- - 1+)' + 2{>^'- + >^l)il- 7+)' + {>^l ^^)'] 

~ 16:7W(A_,-A,)2 (E21) 

Thus we see that both Cg^ are quite small, so that in this case again, not only is the damping of the oscillations 
quite rapid compared to the general decay of the probabilities, but they have a very small amplitude so that the 
oscillations takes place against are very large "background" . 

APPENDIX F: MUTUAL COHERENCE REGIME: OVERDAMPED PLUS UNDERDAMPED 

The Laplace inversion in the Overdamped plus Underdamped case is very similar in spirit to the fully Underdamped 



case. The global denominator Eq. (5.17) is given by 

(A + r2)((2ai + A)2 + J^) + Al{2ai + X) (FX) 



and by defining a = ai — ^2/2, the roots for the variable A + r2 are exactly given by Eq. ( ^.35 ) with a replaced by a 
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In the text, we have argued that we could separate the individual denominator, Eq. ( 5.1^ ) into the individual 
products (A + (72)(A + gi), ie., we neglect the correlation between the two spins when calculating decay rates and 
oscillation frequency. It is then clear that we can directly use Eq. (5.35) for the underdamped spin. 

The decay rates r2 + and r2 + 7 are fully expressed as 



r2 + 7fl; = lR,s + ^2 



iR,s + r2 



and 



r2 + 7 = 2a - i7fl,s(a, J\^l) ^ 2a - (a ~ r2/2)^ 

where the last forms are in the limit J' ^ Ai. The coefficient C2 is given by 
-1 1 



Co = 



(7-R + r2) ((71 -irY + 

which, in the limit ^ Ai reduces to 



4 



r2A+((2ai - 7fl)2 + J^) + -i(8aiA+ - 7^ + r2^/2r) 



c, 



2^+ 



IR,. 



(F2) 



(F3) 



(F4) 



(F5) 



which is of 0(1) is r2 ^ 7fl„s and negligible otherwise. The coefficients relating to the terms oscillating at a frequency 
/X are 



^3 = It^,[NiD2 + N2D^] 
H \D^\ 



and 



1 1 



■[NiDi - N2D2] 



given in a form similar to the coefficients in the underdamped case. The various factors are: 

iVi - r2A+[(2a - 7)2 -fl^ + J'] + ^{8aiA+ - 71 + T2J/2T) ^ ^A'^HaJ/T - T2] 



(F6) 



(F7) 



(F8) 



A2 1 _ 

N2 = -2T2A+ft{2ai - 71) + -^lyi ^ -A^J 



The remaining terms are given by 



^i = (7-r2)(7i?,-7)]+Ai'^:^' 



(F9) 



(FIO) 



D2 = /2(27 - 7fl + r2) ~ 2a^ 

\D'\ - [(7 + ^2? + fi'miR - 71)' + m'] - J' 

The coefficients C5 = C(r2,7i,_R) and Cg = C(7i__R,r2) where the function C{x,y) is defined as 



C{x,y)^ 



1 



r2A+((2ai - x)2 + J^) + ^Alix + 8aiA_ 



{y-x) ((71 - a;)2 + 1/2) 
The coefRcients of the oscillating terms are defined similarly to the Equivalent Spin case 



C7 



1 1 



:[N3D2 + NiDi] 



1 1 



i/i 7f + i^f U'^ vi 



(FU) 
(F12) 

(F13) 
(F14) 
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<^8 = — 2 ! 2 7^[^3£'l - A^4£'2] + --^[A^lDi - 7V2£)2] (F15) 

7i + '^i -D 1^1 
With numerators and denominators given by 

m = [(7i,iJ - lif + i^lW2 - lif + vl] (F16) 

TVi = T2A+[{2ai - 71)^ - + + ^(8aiA+ - 71) (F17) 

A? 

7V2 = -2r2^+z/i(2ai - 71) + -^vx (F18) 

A^3 = - 7i(2ai - 71)] (F19) 

Ni = T2^\viai (F20) 
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Figure Captions 



Figl: The probabilities P^^^-^it) for a system of spins, coupled by Jo, to start in the state | tt) at time t = 0, 
and finish at time t in state \T1T2). We assume A2 = 2.5Ai, and a fairly strong coupling |Jo| = 5Ai. PT^T^(t) is 
plotted as a function of i/Ai. Oscillations between the 2 states | tt) and | [[) occur at the slow "beat" frequency 
A1A2/2IJ0I, with very weak high frequency oscillations superposed, coming from the weak mixing with the states 
I ti) and I It)- These high-frequency oscillations dominate -Pfj,(t) and Pi'\{t), albeit with amplitude reduced by a 
factor ~ (Ai + A2)/2Jo. 



Fig2: The same probabilities as in Fig. 2, but now for weak coupling; we have A2 = 2.5Ai again, but now 
I Jol = Ai/10. 



Fig3: Diagrammatic interpretation of the equation of motion for the density matrix. We show equation 4.2 in 
Feynman diagram form; J is the propagator (or Green function) for p. 



Fig4: Contributions to the single spin influence functional. The label a = 1,2 labels one of the 2 spins, and 
7q, = Fq, — j$Q;, (see Eq. (4.15) and (4.16)) These contributions are exponentiated to give the influence functional Eq. 
(4.22). 



Fig5: Contributions to the interaction functional F12 (Eq. 4.23)). There are only interactions between the path 
associated with the spin 1 and that of the spin 2. Here, 712 = ri2 — i$i2, as given by Eq. (4.17) and (4.18). 



Fig6: Example of a typical path for the 2-spin system. The top path is for the 1st spin, and shows transitions 
between eigenstates of rf , at times ty, the bottom path is for the 2nd spin, with transitions at times Uk- 



Fig7: The various regimes in which the PISCES model dynamics are solved for analytically in the text. We assume 

in the figure that the 2 spins are equivalent, ie., that A* ~ A^ = A2, and ai = a2 = a. We also assume weak 
damping, so that a <C 1. If a is much larger, ie., a ~ 0{1), there is no mutual coherence phase, and motion in the 
locked phase (where the 2 spins rotate rigidly together) is overdamped at any T. We do not discuss the perturbative 
regime in this paper - in this regime correlations are very weak between the 2 spins. If the 2 spins are different (ie., 
Oil 7^ ct2, and/or AJ 7^ A2), then one must draw 2 separates diagrams of this type, one for each spin - then it is 
possible for the 2 spins to be in different phases. The Mutual Coherence regime is strictly defined by the conditions 
A* /a '3> T ':§> {j"^ + A*^)^/^, but since we always assume ^ A* in this paper (ie.. well away from the perturbative 
regime), the inequality T ^ {j"^ + A*^)^/^ is equivalent to the condition T ^ J used in the text. 



Fig8: The various dynamical regimes which describe the motion of the PISCES system in the locked phase, as 
functions of T and of the coupling ac between the locked 2-spin complex and the sea of oscillators. The overdamped 
relaxation phase is divided into 2 parts by the "Toulouse line" a,, = 1/2. On the "weak relaxation" side, the damping 
rate Fg decreases with increasing T, and is small; on the "strong relaxation" side, is large and increases with 
increasing T. When T = 0, the motion is still damped for any finite ac (with oscillations for ac < 1/2), but for 
ac > Ij the system is completely frozen by the coupling to the oscillators. 



Fig9: Probabilities of occupation for a system in the "correlated relaxation" regime, with parameters ai = 1.5, 
a2 = 2, at a temperature such that T/Ai = 100, T/A2 = 300 and a ferromagnetic coupling J /2T = —0.02; the time 
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is expressed in units of inverse temperature. The system starts at i = in the state | tt); Paf3 is then the probabiHty 
that at time t the system has rf = a and t| = f3. We plot the probabihties as a function of ln{t) in order to clearly 
show the different relaxation times involved. 



FiglO: Imaginary part of the Fourier transform of P]i {t) in the mutual coherence case. The peaks represent 
the oscillation frequencies, and their width is proportional to their damping rates. We use the values A/ J' = 0.6, 
J'/T = 0.1, T = 100 (in arbitrary units) and a = iraT = 0.25. The real frequency iX is in the same units as the 
temperature. 



Figll: Imaginary part of the Fourier transform of P]'i{t) in the Overdamped plus Underdamped case. The param- 
eters are Ai/J = 0.6, A2/J = 0.1, ai = 0.25 and a2 = 15 (where = napT). Again J /T = 0.1, T = 100. Notice 
that we still have 0:2 <C 1, although the spin 2 is overdamped. The strong broadening of the oscillations caused by 
the environment does not allow a clear separation of the two different oscillation frequencies. 



Figl2: Paths for P|| with ni = 7^2 = 1 (cf. equation (A2); there are 6 possible orderings of the transitions, which 
must be summed over. In Fig 12(a) we also show the states associated with the periods between the transitions. 



Figl3: Schematic representation of the path integral for P||. Wc sum over all possible chains of coupled clusters. 
The charges of the boundary sojourns are fixed at = +1 by boundary conditions and there is a summation over the 
charges of all the internal sojourns. 



Figl4: Representation of a cluster, used for the summation of Ptit2- this graph, the charge of the lower path 
sojourn is fixed and constant. The charges of the boundary sojourns of the upper pat hs a re fixed, and a summation 
is done over the charges of the internal sojourns. This summation then results in Eq. ( |D5| ). 
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